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ABSTRACT 


Steam condensation heat transfer measurements were made 
in a 5-tube test condenser having an additional perforated 
tube to simulate up to 30 active tubes. Results were 
obtained for smooth tubes and roped tubes wrapped with wire. 
A Sieder-Tate equation was used to correlate the inside 
heat-transfer coefficient. For smooth tubes, a leading 
coefficient of 0.029 was found, while it was 0.061 for the 
roped tubes. The average condensing coefficient measured 
for 30 smooth tubes was 0.59 times the Nusselt coefficient 
calculated for the first tube. When the smooth tubes were 
wrapped with wire, this ratio increased up to 0.86. 
Further, roped tubes without wire experienced a ratio of 
0.63, while roped tubes wrapped with wire resulted in a 
ratio of 0.86. These preliminary data show that wire- 
wrapped tubes may lead to a significant reduction in 


condenser surface area. 
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1. AL BA ROUN 


In recent years, there has been a continued interest in 
the reduction of the size and weight of propulsion systems 
aboard both surface vessels and submarines. Especially, 
actual dimensions of naval condensers have a critical bear- 
ing on cost and performance of the ship. Often, compactness 
1s more important than thermal effectiveness when the over- 
all performance of the ship is considered; in a submarine, 
the diameter of the pressure hull can depend on the dimen- 
sions of the condenser. 

The importance of compactness justifies measures to 
raise the overall heat-transfer coefficient of condenser 
tubes despite the penalties which may occur, i.e., the 
increased pumping power and tube cost. 

Naval condenser design is based upon the Heat Exchange 
Institute (HEI) specifications for steam condensers [Ref. 1] 
and also the standards of the Tubular Exchange Manufacturers 
Association (TEMA) [Ref. 2]. Search [Ref. 3] investigated 
the present condenser design processes, including the 
feasibility of enhanced heat transfer in naval condensers. 
He concluded that the current design is very conservative, 
and he predicted that a forty-percent reduction in condenser 
weight and volume could be achieved depending on the heat- 


transfer enhancement method used. 
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In recent years, many research efforts have been 
directed to the study of heat-transfer enhancement tech- 
niques and their application to heat-exchanger design. Webb 
[Ref. 4] has summarized extensive works of augmentation 
techniques. At the Naval Postgraduate School, Beck [Ref. 
5], Pence [Ref. 6], Reilly [Ref. 7], Fenner [Ref. 8] and 
Ciftci [Ref. 9] conducted experimental research into various 
kinds of enhancement schemes employing a single-tube test 
condenser. 

The above-mentioned investigations concluded that, for 
the same diameter tube, the overall heat-transfer coeffi- 
cient of enhanced tubes can exceed those for smooth tubes by 
almost 100 percent. Reilly and Fenner [Ref. 19] revealed 
that most of the above mentioned augmentation occurred on 
the cooling-water side due to a combination of increased 
Surface area, and increased turbulence and swirl in the 
cooling water flow. Little or no improvement occurred on 
the steam side. Eissenberg [Ref. 11] performed an extensive 
study on condenser-tube, heat-transfer coefficients uSing a 
multi-tube bundle. 

In order to investigate the outside heat-transfer 
performance of various enhanced tubes in tube bundles, 
research was conducted at the Naval Postgraduate School. 
Noftz [Ref. 12] modified a test apparatus initially designed 


by Morrison [Ref. 13] to simulate a tube bundle using five 
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active tubes arranged in a vertical plane. A perforated 
tube was located at the top of the bundle, through which 
water was flooded to simulate bundles having up to 38 tubes 
in a vertical row. His investigations determined that the 
heat-transfer coefficients for a given tube of the tube 
bundle increased as the mean vapor velocity increased, but 
decreased as the amount of condensate inundation increased. 
The experimentally found values for the heat-transfer 
coefficients were comparative with the Nusselt theory. 
Based upon research done currently, it is evident that 
present day smooth-tube steam condensers, operating under 
typical conditions, have limitations in their thermal 
efficiency, due to a large thermal resistance which occurs 
on the tube side of the condenser. This resistance is 
generally larger than any of the thermal resistances that 
occur on the steam side, in the tube wall, those due to 
fouling, or due to noncondensable gases. However, employing 
enhanced tubes, the inside heat-transfer coefficient can be 
increased by 100 to 200 percent over the smooth-tube case. 
The outside heat-transfer coefficient, on the other hand, is 
increased by only 10 to 50 percent. Іп this situation, the 
thermal resistances on the inside and outside of the tube 


can be approximately equal. 
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Webb [Ref. 4] reported that the dominant thermal resis- 
tance in film condensation is that of conduction across the 
condensate film and, therefore, a surface geometry that 
promotes reduced film thickness will provide enhancement. 

Thomas, et al, [Ref. 14] tested ammonia condensation on 
a smooth tube with a wire wrapped in a helical manner. The 
measured condensing coefficient was approximately three 
times that predicted by the Nusselt equation for a smooth 
tube.  Surface-tension forces draw the condensate to the 
base of the wires, which act as condensate run off channels. 

Webb [Ref. 4] stated that, when noncondensables are 
present, an additional thermal resistance is introduced in 
the gas at the vapor-liquid interface. Mixing in the gas 
film will substantially reduce this thermal resistance. 
Therefore, the maintenance of high vapor velocities, or 
Special surface geometries that promote a higher heat- 
transfer coefficient in the gas film will substantially 
alleviate the performance deterioration due to 
noncondensables. 

Cunningham [Ref. 15] presented in his paper that, for 
the roped tubes on the vapor (shell) side, the enhancement 
is achieved by improved condensate drainage, while on the 
coolant (tube) side the helical ridges increased turbulence 


and, as a result, the inside convective coefficient. 
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Improvements on the condensing side up to 199 percent 


have been reported for single-tube tests [Refs. 16;17]. 


Although titanium has a low thermal conductivity, it pro- 


vides a high resistance to erosion and water-side fouling. 


Titanium tubes with enhancement both inside and outside are 


commercially available through the Wolverine Tube Division 


of Universal Oil Products, Inc. The applications of these 


tubes having all the inherited properties of titanium are 


promising for naval condensers. 


The goals of this thesis were therefore to: 


l. 


Obtain baseline heat-transfer performance data for 
the test condenser utilizing 16 mm O.D. smooth 
titanium tubes. 

Conduct steam condensation tests with the following 
enhanced tube geometries to determine steam-side 
heat-transfer coefficients in relation to smooth- 
tube performance: 

a. Wolverine "roped" tubes. 

b. Wolverine "roped" wrapped with titanium wire. 


c. Smooth tubes wrapped with titanium wire. 
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II. THEORETICAL BACKGROUND 


The combined effect of vapor shear and inundation on the 
condensate film heat-transfer coefficient for cylindrical, 
horizontal tubes within tube bundles is a very complex and 
still insufficiently-understood subject, which is of impor- 
tance to the efficient design of steam condensers. Although 
many researchers have studied this subject both theoretic- 
ally and experimentally, there is no accurate me thodology 
available for predicting the condensate-film, heat-transfer 
coefficient within tube bundles. 

In 1916, Nusselt conducted his pioneering analysis for 
the simple case of condensation occurring on the outside of 
a single, isolated, horizontal tube. He idealized the 
problem by making the following assumptions for single tubes 
as stated by Nobbs [Ref. 18]. 

1. The wall temperature is constant. 

2. The flow is laminar in the condensate film. 


3. Heat transfer in the condensate is by conduction, 
and subcooling may be neglected. 


4. The fluid properties are constant within the 
condensate film. 


5. The forces due to hydrostatic pressure, surface 
tension, inertia, and vapor-liquid interfacial 
shear are negligible when compared to the viscous 
and gravitational forces. 


6. The surrounding steam and vapor/liquid interface 
are at saturation temperature. 
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7. The film thickness is small when compared with 
normal tube diameters and the effects of curvature 
are small. 


Based on the above assumptions, Nusselt predicted the famous 
relationship for the heat-transfer coefficient: 


1/4 


3 
k p- ° 
h ЕО 725 T 
Nu j UD (T -T ) 
sat `w 


(21) 


In order to simulate and analyze a tube bundle, 
Eissenberg [Ref. 11] stated the following additional 
assumptions: 

8. Condensate drains as a laminar sheet from a tube on 
to the tube directly underneath in such a way that 
velocity and temperature gradients are not lost in 
the fall between tubes. 


9, The saturation temperature and the tube-wall 
temperature are constant for all tubes in the bank. 


Jakob [Ref. 19] extended the Nusselt analysis for 
filmwise condensation heat transfer on a vertical in-line 
row of horizontal tubes as shown in Figure la. 

The above-mentioned assumptions were combined with the 
assumption of constant temperature drop across the 
condensate film for all the tubes, and the average 
coefficient for a vertical row of N tubes was predicted to 


be: 


1/4 


_ k?o(p-P,)h. ‘а 
ЕҢ ааа пъ 
Ay ен 25 UN D (T т) 
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Figure la. 


Idealized Condensation 
on Banks of Tubes 





Fiqure Ip- 


More Realistic Picture of 
Condensation on Banks of Tubes 
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Upon, dividing equation (2.2) by equation (2.1), the Nusselt 
theory can be expressed as: 


ћ = Деда 
Ny / (2.3) 
Nu 


Equation (2.3) can be also expressed in terms of the 


local coefficient for the N-th tube: 


(N-1) 3/4 (2.4) 

In reality, condensate does not drop off in a continuous 
laminar sheet, but drops off instead by discrete droplets of 
liquid, as shown in Figure lb, depending upon the surface 
tension of the condensate. These droplets create ripples in 
the condensate film, and thereby decrease the performance 
degradation due to inundation. Based on his research, Kern 
[Ref. 20] proposed a less conservative relationship: 


-1/6 


a =N (255) 


hs 


Or, in terms of the local coefficient for the N-th tube, 


Sue mery W SŠ (2.6) 


23 





Chen [Ref. 21] considered the following conditions: 


1. the momentum gain of the falling condensate between 
tubes, and 


2. the condensation of vapor on the condensate between 
tubes, 


and concluded that: 


E 1/4 
ШІ -1/4 E 1-0.685 + 0.025€ 
Nu 
where; 
CHAT 
Еш = / C= —— 
fg fg 
and EN 
E. 2. 
5 P 
Е 


The above approximate expression, due to Chen, is valid for 
most ordinary applications. 

Experimental work doen by Eissenberg [Ref. 11], in order 
to investigate the effects of steam velocity, condensate 
inundation. and noncondensable gases on the heat-transfer 
coefficient, revealed that condensate does not always drain 
onto tubes aligned vertically, but can be diverted sideways, 
caused by local, vapor-flow conditions. The condensate thus 
follows a staqgered path as shown in Figure 2. 

Eissenberg, making the assumption that the flow is 
dominated by gravity, stated in this side-drainage model, 


that condensate strikes the lower tubes on their sides 
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rather than their tops. Therefore, the inundation effects 
influence the condensate flow only on the lower half of the 
tubes, which transfer less heat than the upper half. Based 
on the above-mentioned model, Eissenberg obtained the 
following formula: 


B. -1/4 


—— = 0.60 + 0.42 N (2.8) 


"Nu 


Much experimental research has been conducted for 
studying the effect of condensate inundation. In general, 
the obtained data are highly scattered as shown in Figure 3. 
Berman [Ref. 22] made a comprehensive comparison of filmwise 
condensation data on bundles of horizontal tubes, and he 
concluded that the wide variation in expermental data for 
tube bundle inundation is caused by the following variables: 

l. bundle geometry (in-line or staggered), 

2. tube SESTO, 

3. type of condensing fluid, 

4. operating pressure, 

5. heat flux, and 

6. local vapor velocity. 

In addition. noncondensable gases, and insufficient steam 


for lower tubes can cause such scattering of data. 
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Nobbs [Ref. 18] used one active tube in a dummy tube 


bundle. 


He simulated additional condensate by using three 


porous tubes. Based on his results, he concluded the 


following: 


b 


Vapor velocity increases the condensate heat- 
transfer coefficient on both inundated and 
uninundated tubes in a given tube bundle. 


The effect of inundation is to reduce the heat- 
transfer coefficient. 


The condensate drainage path is often not vertically 
downwards but in a diagonal direction. This can 
result in tubes receiving different amounts of 
inundation. 


Marto and Nunn [Ref. 23] made a comprehensive survey of 


the effects of noncondensable gases. They stated that these 


can be classified into one of two categories: 


e 


the introduction of an additional local thermal 
resistance due to the propagation of noncondensable 
gases towards a condensing tube surface under the 
influence of a gradient in partial pressure, and 


the cumulative effect of gas blanketing where uneven 
rates of condensation in a condenser bundle eventu- 
ally lead to regions where tubes are inoperative in 
a condensing role. 


Experiments done recently by Cunningham [Ref. 15] revealed 


that noncondensable gases have a less effect on the perform- 


ance of finned tubes compared to smooth tubes. Based on the 


above discussion, it is clear that the performance degrada- 


tion due to noncondensables must be taken into consideration 


in realistic designs of tube bundles. 
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Another very important factor in the performance of a 
condenser is the effect of vapor velocity. As noted 
earlier, the vapor shear plays a beneficial role. Berman 
and Tumanov [Ref. 24] conducted experiments on a single 
zontal tube placed in a bank of uncooled neighboring 
tubes. For vapor in vertical downflow, they found a 


relation between the vapor Reynolds number, and the heat 


flux as: 
h = ШЕВАИ 
= = 1+ 9.5 x 107? (Rev) 
N 
with the restriction, that ju Mn 
(Rev) < 50 


Eissenberg [Ref. 25] has stated that in designing exper- 
imental bundles, the combined effects of inundation and 
vapor shear are very important. In fact, the use of narrow 
condenser bundles to experimentally study vapor shear and 
inundation effects is preferred to simulate large conden- 
sers. However, small narrow tube bundles can create errors, 
due to the following causes: 


l. Wall flow: condensate drainage may reach side 
walls; 


2. Dummy tubes: condensate inundation may disperse; 


3. Vapor lanes: vapor may bypass along the side walls; 
and 


4. Noncondensable gases: they can affect condensation 


even at low concentration in the bulk stream, 
particularly if steam is recycled. 


29 








III. EXPERIMENTA А Т 


E TEST FACILITY 

The test facility shown in Figures 4 and 5, was designed 
and built by Morrison [Ref. 13] and modified by Noftz 
[Ref. 12] to simulate an active tube column having up to 30 
tubes in increments of five tubes deep (i.e., five, ten, 
fifteen, etc.). Some elements of the original test facility 
were modified to allow for more efficient operation of the 
facility. 

A detailed description of the components used in the 
test facility is given in Reference 12. Only a short 
description of these components will be found in this 
report. Particular attention. however, will be focused on 
the experimental tubes. Calibration procedures for 
components requiring calibration are outlined by Reilly 


[Ref. 7]. 


B. STEAM SYSTEM 

The steam system shown in Figure 6 was modified from 
Noftz's initial design. From the house supply line, steam 
flows via a 19 mm O.D. stainless steel line through a steam 
supply valve (MS-3) to a cast iron steam separator. The 
steam continues through the system past two Nupro bellows 


valves, which were used in conjunction with the supply valve 
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to regulate the steam supply pressure. From these valves, 
the steam flows into the test condenser diffuser. 

The steam supply pressure was monitored by a pressure 
gage just downstream of the steam supply valve; also a 
compound gage just after the Nupro valves was used to 
monitor the pressure in this line. 

The operator had no control over the state point, 
quality or noncondensable gas content since house steam was 
used. However, especially at nights and during weekends, 
the state point of the steam at the inlet to the test 


condenser was found to be nearly constant. 


C. TEST CONDENSER 

The test condenser shown in Figure 7 was unchanged from 
Noftz's initial design. Steam enters via the top, passes 
through the transition piece, the vortex annihilator, and 
finally through the diffuser. The dimensions of the test 
condenser were 305 mm X 305 mm X 79 mm and it was made of 
stainless steel. These dimensions allowed for a maximum of 
twenty-seven l6 mm O.D. tubes arranged in an in-line config- 
uration of three columns of nine tubes each. For this 
experiment, the in-line configuration was used with a middle 
column of five active tubes flanked on either side by a 
column of five dummy tubes. Just on top of the upper active 
tube, a perforated, distilled water supply tube was posi- 


tioned, and flanked on both sides by dummy tubes. In order 
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to conduct the experiments, a square, in-line arrangement of 
tubes, with a pitch-to-diameter ratio of 1.5 was used. 

A vertical slot (Figure 7) along the centerline of each 
condenser end plate was used for active and perforated tube 
installation. The tubes were positioned using nylon tube 
sheets that were attached to the exterior of the condenser 
end plates. 

To minimize heat losses, and also to prevent leaks from 
the tube sheets, each test condenser side was provided with 
a nylon tube sheet one-inch thick with six holes 
(S/D = 1.5), having about 0.5 mm tube clearance, which 
allowed the tubes to be easily slid into the test condenser 
through the tube sheets. The exterior side of each tube 
sheet had grooves to support O-rings. The aluminum tube 
sheets had six holes, having a tube clearance of about 
0.5 mm and were used as sealing plates. The diffuser, 
exhaust plenum, transition piece. vortex annihilator, and 
the exhaust piping which are shown in Figures 8 and 9, were 
insulated with rubber insulation. 

A viewing window allowed viewing of the condensation 
process. A double-walled glass window was used, and heated 
air was fed through the clearance between the two glasses to 


eliminate fogging on the inner glass. 
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D. TEST CONDENSER TUBES 

Four kinds of active tubes were tested in this experi- 
ment. The tubes were manufactured by the Wolverine Division 
of Universal Oil Products, and all were made of titanium. 

The first kind of tubes, shown in Figure 198 were smooth 
titanium tubes of 16 mm O.D. with a 1.65 mm wall thickness. 
The second kind of tubes, shown in Figure ll, were single- 
start, helically-corrugated tubes, designated as Low 
Pressure Drop (LPD), with 16 mm O.D. and a 1.65 mm wall 
thickness. The third kind of tubes, shown in Figure 12, 
were also single-start, helically-corrigated LPD tubes. 
wrapped with titanium wire of 1.58 mm O.D. The fourth kind 
of tubes, shown in Figure 13, were smooth tubes, wrapped 
with 1.58 mm O.D. titanium wire on the same pitch as the 
roped Wolverine tubes. In order to wrap the titanium wire, 
the tubes were attached into a lathe, and the wire was 
welded at the start of the helical groove. Then, under a 
constant tension of 4 kg, the wire was wrapped manually by 
turning the chuck of the lathe. Finally, the free end of 
the titanium wire was welded at the end of the helical 
groove. 

The 660-mm-long active tubes were connected to separate 
cooling water supply and discharge lines. The dummy tubes 
Flanking the active tubes were made of 16 mm O.D. stainless 


steel. These tubes served to direct the steam flow so as to 
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Photograph of Wire Wrapped Smooth Titanium Tube 


Figure 13. 





Simulate actual conditions in a condenser. Cooling water 
was not supplied to these tubes and they did not penetrate 
the test condenser end plates. Special characteristics of 


the Wolverine tubes are listed in Table I. 


E.  PERFORATED TUBE 

The perforated tube water supply system is shown in 
Figure 14. This system consisted of a perforated tube 
(located above the uppermost active tube), a water heater 
which served as a supply tank, a rotameter in order to 
control the amoung of water supplied to the perforated tube, 
a pump driven by a 1/2-HP electric motor, a condensate pump 
and associated piping system and valves. The active length 
of the copper-nickel, perforated tube was 305 mm, which was 
identical to the length of the test condenser. Supply water 
entered one end of this tube; the other end was sealed off. 

The supply tank was unchanged from Noftz's initial 
design, except for the addition of a thermocouple in order 
to have a direct indication of the exact temperature of the 
water supplied to the perforated tube. The water, heated to 
the temperature of the condensate leaving the bottom tube, 
was fed to the perforated tube by a 1/2-HP, electric-motor- 
driven pump, via a rotameter and a recirculation valve. The 
flow rate to the perforated tube was controlled by using the 
rotometer and the valve provided in the heater-water 


recirculation line. Also, another modification was made to 


44 





TABLE I 


DESIGN CHARACTERISTICS OF WOLVERINE 
CORODENSE TUBES (TYPE LPD) 


Outside diameter (in) 2 4573 
Wall thickness (in) 779 7102 
Outside diameter (ft) : 0.0520 


Heat Transfer Surface and Ratio 


a. Outside A. (ft^/£t) 0163 
b. Outside to inside A/A, : 1.141 
Inside diameter D; (ft) : 0.0435 
Cross section for flow Ans (£t^) : 0.00163 
Number of groove starts monec) 
Pitch (in) 22505 00 
Groove radius (in) : 1/32 nominal 


Depth-transition length 
(next to plain section) (in) : 1-3/4 max 
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feed the condensate collected in the hotwell, via the 
condensate pump, to the supply tank (heater) in order to 
facilitate the inundation of up to 30 tubes under atmos- 
pheric conditions. 

The main feature of the perforated tube was to inundate 
with condensate from above, in order to simulate a tube 
column of more than five active tubes. When the wrapped 
Wolverine tubes and the smooth wrapped tubes were tested, 
the perforated tube was also wrapped with titanium wire, by 


using the above-mentioned technique. 


F. CONDENSATE SYSTEM 

The condensate system shown in Figure 15 was modified 
from Noftz's initial design, and consisted of the test 
condenser and hotwell, the condensate pump, piping and 
valves. 

The test condenser hotwell collected the condensate 
produced by the test tubes. With valve C-1 closed, the 
condensate mass flow rate from the test condenser could be 
measured using a stop watch. Opening the valve, the 
condensate was fed, via the condensate pump, to the supply 
tank for the perforated tube, or dumped into the building's 
drainage system. The condensate line connecting the test 
condenser and the test condenser hotwell were insulated 


using rubber insulation. 
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С.  COOLING-WATER SYSTEM 

The cooling-water system was a partially-closed system 
as shown in Figure 14, and it was unchanged from Noftz's 
initial deSign. House water was used for the test facility. 
Cooling water was stored in a 1.2 meter, cubical, plexiglass 
supply tank, and was pumped by a 5-HP, electrically-driven 
pump via 51-mm-O.D. plastic piping to a manifold. Five 
rotameters were attached to the manifold to measure the flow 
rate through each active tube. Five regulating valves were 
used to obtain any water velocity between @ and 5 m/s within 
the tubes. The cooling water passing through the rotameters 
was fed via 16 mm O.D. stainless steel tubing, and tygon 
tubing to the active tubes. The total tube lengths were 
long enough (over 2 meters) to ensure hydraulically 
fully-developed flow, and no swirling into the tubes. After 
flowing through the tubes, the cooling water was passed 
through mixing chambers, where the temperature profile was 
destroyed, to facilitate the steady measurement of the 
outlet cooling water temperatures. After the mixing 
chambers, the cooling water was collected in the supply 
tank. A 7.5-HP, electrically-driven pump was used to pump 
water from the supply tank, through a filter, to a cooling 
tower in order to minimize temperature rise at the water 


inlet. 
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The cooling tower was located outside the building and 
was composed of four truck radiators across which air was 
blown by means of a fan. The entire system, consisting of 
the heat exchanger and the fan was enclosed in a wooden 
structure with louvered openings to provide enough 
ventilation. 

The tygon tubing was secured to the inlet and outlet 
sections of the active tubes by means of hose clamps. The 
outlet sections, including the mixing chambers, were 


insulated with rubber insulation. 


Н.  INSTRUMENTATION 
1. Flow Rates 

a. Foulton rotameters were used to measure the flow 
rate of cooling water for each active tube. Starting from 
the top active tube to the bottom one, the rotameters were 
calibrated giving 199% maximum flow rates of 66.941, 
72.641, 71.541, 72.8+1, and 73.341 kg/min. 

b. The perforated tube water rotameter was 
calibrated giving a 100 percent volumetric flow rate of 
9694-30 ml/min. 

C. All rotameters were calibrated using the 


procedures noted in Appendix A of Reference 7. 
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2. Temperature 

Stainless-steel sheathed, copper-constantan 
thermocouples were used as the primary temperature 
monitoring devices. Three thermocouples were used to 
measure the inlet cooling water temperature; and twelve 
thermocouples were utilized to measure the outlet cooling 
water temperature. 

Additionally, two thermocouples were used to measure 
the steam saturation temperature; one thermocouple was used 
to monitor the condensate temperature, and one was used to 
measure the vapor temperature. Table II lists the locations 
monitored. 

A Gulton Industries, West 20, 0-5009 Е temper- 
ature controller was used to regulate the temperature of the 
perforated tube supply water. The controller had a 
manufacturer-stated accuracy of %0.5% of the span (about 
+1.259 Е). It was impossible to obtain the precise 
temperature for the perforated tube water, because of the 
low response of the heating system. 

During inundation runs, the condensate temperature 
rose with the increasing number of tubes. For example, the 
condensate temperature after the 5th tube was 94.69 C 


and it increased to 97.65? C after the 30th tube. 
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3. Pressure 

Several different types of pressure measurement 
devices were used in this facility. They were: a Bourdon 
tube pressure gage which was used to measure the steam 
supply pressure, located downstream of the steam supply 
valve; an absolute pressure transducer which was used to 
measure the test condenser pressure; and, a compound 
pressure gage was also used to measure the pressure drop 
downstream of the Nupro valves. The pressure transducer was 
calibrated against a mercury manometer. 

4. Data Collection and Display 

A Hewlett-Packard HP-3054A Automatic Data 
Acquisition System, with a HP-9826 computer and a HP-2671G 
printer were used to record and display the thermocouple and 
pressure transducer readings. The temperatures were 
recorded in degrees Celsius, while the pressure was recorded 
in volts which was then converted, using the calibration 
curve, into mm Hg absolute. The pressure transducer was 
assigned to channel 019 of the HP Automatic Data Acquisition 
System, while the thermocouples were assigned channels as 


indicated in Table II. 
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TABLE II 


CHANNEL NUMBERS FOR COPPER-CONSTANTAN THERMOCOUPLES 


Location Channel 
Toi #1 000 
T #3 001 
Tei #5 002 
ед #1 003 
Та #1 004 
Tas #1 005 
т #2 006 

со 
- 12 007 
ТЕЬ # З 008 
па. # З 009 
E #4 010 
co 
Tag #4 011 
T #4 012 
CO 
Dass #5 013 
Das #5 014 
Tessa: 015 
a 016 
Tod 017 
р 018 
уар 


DE 
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IV. PROCEDURES 


A. INSTALLATION AND OPERATING PROCEDURES 


ШЕ. 


Preparation of Condenser Tubes 


Prior to installation. the titanium tubes were 


cleaned using a chemical cleaning method [Ref. 26]. The 


steps in this cleaning process are as follows: 


d. 


b. 


e. 


f: 


Swab the tube surface with acetone to remove grease. 


Using a test tube brush, brush the inside surface of 
the tube with a 50$ sulfuric acid solution in order 
to remove any oxides. Also. apply this solution 

to the outside surface of the tube. 


Rinse the inside and outside of the tube with tap 
water. 


Apply. using a brush. a 50$ solution of sodium 
hydroxide mixed with an equal amount of ethyl 
alcohol. at the boiling temperature (about 

859 C) to the outside surface of the tube. 
Rinse the tube with tap water. 


Rinse thoroughly with distilled water. 


Prior to any run. the condenser tubes had to be prepared to 


ensure filmwise condensation- Exterior and interior 


surfaces were cleaned to ensure proper wetting character- 


istics. 


Also, the tubes were cleaned by running steam at 


atmospheric pressure through the test condenser for about 


twenty minutes without the cooling water running through the 


tubes. 
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It was also found that, when drop-wise condensation 
occurred, rinsing the tubes using the perforated-tube water 
was sufficient to restore film-wise condensation. 

2. System Operation and Steady-State Conditions 

Complete operating instructions are listed in 
Appendix A. The steady-state condition was reached about 
three hours after initial system light-off, and about 
fifteen minutes after changes to the cooling water or 
perforated tube supply water flow. 

When a steady-state condition was reached. the runs 
were made. The duration of each run was approximately one 
minute. For each condition. for example, for tubes 11 
through 15. five consecutive runs were made and average 
values were computed. 

The following data were taken automatically by the 
data-acquisition system: 

a. the thermocouple readings, апа 

b. the pressure transducer reading. 
Also. the following data were read into the computer, 
through the keyboard, for each run: 

а. the setting of each rotameter, and 


b. the test condenser hotwell levels. 
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B. DATA REDUCTION PROCEDURES 


l. Overall Heat-Transfer Coefficient 


The heat-transfer rate to the cooling water is given 


О = НС wy Тео" Тел) | (451) 


The heat-transfer rate can also be found from the 


overall heat-transfer coefficient by: 
Q-U A * LMTD (ДҮР 2!) 
O O 


where 


Tear” Taj? "5 (T 57 To 9! 


Тс 
(ване) (4.3) 


sat СО 


LMTD - 


After combining equations (4.1). (4.2) and (4.3), it 
is found that 


С W — 
ee 93. 


Е О sat со (4.4) 


2. inside Heat-Transfer Coefficient 
The heat-transfer coefficient on the inside is 
calculated from the Sieder-Tate relationship as described in 


Holman [Ref. 27]: 
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In the above equation. с. is referred to as the 
Sieder-Tate coefficient. The remainder of the right-hand 


side of the above equation. 


is referred to as the Sieder-Tate parameter, X. 

For thermally and hydrodynamically developed flow in 
tubes, the Sieder-Tate coefficient equals 0.027. However, 
when short tubes are considered, as in the case of this 
experiment, fully-developed conditions are not attained. 
Therefore, the value of Ci must be found experimentally. 

The Wilson plot was used to arrive at the value of 
the Sieder-Tate coefficient. The Wilson plot was developed 
in 1915 [Ref. 281. It is merely a plot of 1/0, versus the 
inverse of the Sieder-Tate parameter (which is proportional 
to the inverse of cooling water velocity raised to the 0.8 
power). The reasoning behind the Wilson plot is shown in 
the following development. 

Consider the equation for the overall heat-transfer 


coefficient: 
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I O l 
=— = — + R + (4.6) 
U A.h; w h. 
Or 
да D. 
7 ~ оте ћ (457) 
UN А. СКХ w " 
where 
MC 
D, in(g-) 
R = = 
Мм 2k 


For the Wilson-plot method to be successful, no 
must be kept constant. When steam condenses on the shell 
side, the above condition can be achieved only if the heat 
flux (q") is kept constant. 

As required for the Wilson plot. when the cooling 
water velocity increases. ni increases. U increases, 
and finally q" increases. 

Consider 

q" = U LMTD (4.8) 

In order to keep q" constant, one must therefore 
decrease LMTD by lowering the steam saturation temperature. 
This requires trial-and-error setting of condenser pressure 


which is a very difficult task. 
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To avoid this difficulty and yet arrive at a Sieder- 


Tate constant not affected by the varying q", a modified 


Wilson-plot method was used. 


This method was developed by 


Wanniarachchi [Ref. 291, and the required steps are listed 


below: 


5А. 


Assume Ci = 0.027 


Calculate 


Dm Tuc 
СО СЪ 


an ( тезі еі) 
Feat “co 


Calculate 


LMTD = 


Q = me, Ne > lei 


Calculate 


U A OG 
О A,° LMTD 


For the first data point only. do the following: 


a. Assume 


X = R 0.8 p 173 ЕШ, 0.14 
= ЈЕ u 
М 
с. Calculate 
uk 0.8 1/3 u 0.14 
h; D. Ci В. P. (m) 
М 


DS 





d. Calculate 


f. Repeat steps a through e until C 
assumed in step a approximately Equals 
Ce calculated in step e. 


а. Calculate 


Fi OS. 
h ЧАШ h Aa. W (4.9) 
O O 143991 
h.  Assign 0, - 0 


NOTE: The second subscript of h on the left-hand side of 
equation (4.9) refers to the first data point. 


5B. For all data points except for first data 
point. do the following: 


a. Calculate 











dl = о-в m WE ^. 
h. Б a 
L Чо ы ШС, a. 
where 
Do UN on Q1 


b. Calculate 
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с. Calculate the Sieder-Tate parameter 


6. Repeat steps 2 through 5 for all data points. 


7. Plot — versus 


l 
А X 
NOTE: It is more reasonable to plot h. versus X. 


However, l/h. versus l/X was plotted in order to 
be consistent with the original Wilson-plot method 


1 M. 


| 
a" versus x 


О 


8. Obtain the slope (m). of the least-squares-fit. 
straight line. 


9. Calculate the Sieder-Tate constant. 





10. Repeat steps 2 through 9 until the assumed and 
calculated C; values are approximately equal. 


A listing of the computer program used in this 
method can be found in Appendix D. 
3. Qutside Heat-Transfer Coefficient 
The heat-transfer coefficient on the outside is 
calculated using the following steps: 
l. Calculate the average bulk water temperature. 


ie = (T 
C 


b + T.,)/2 


О 
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Evaluate the thermophysical properties based on 
the average bulk water temperature. 


Calculate the cooling water velocity. 


Calculate the water-side Reynolds number. 


DE 
A= LS 
w 


Calculate the heat transferred to the cooling 
water. 


pom (Тео Б Пе Cow 


Calculate the heat flux. 


" Q 
а =, 
T Do L 


Calculate the Nusselt coefficient using the 
formula: 


T. head 1/3 
h = 0.651 


T (4.10) 
Nu He, D. d 


NOTE: The above formula was employed since no direct 
measurement of tube wall temperature was made. 


To compare the condensate film temperature, as 


required for equation (4.9), an iterative scheme was used as 


outlined below: 
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a. Assume Te = Tat’ 
b. Evaluate the relevant thermophysical properties. 
which are included in equation (4.142). 


с. Calculate the Nusselt coefficient using 
equation (4.10). 


d. Evaluate Te © using the formula: 


= EN Ak. 
E Tsat Ayu 


Е с web m 


e. Repeat steps a through e until T, assumed 
in step a approximately equals Т; c Caleu- 
lated in step e. = 


8. Calculate the inside heat-transfer coefficient using 
the formula: 


W 0.8 
pitu pa; e Pr = „Се (4.11) 


NOTE: In order to determine the C., the average inner 
wall temperature must be known. As noted earlier, 
this temperature is not known directly and it must be 
found iteratively, as described below: 


a. Assume a correction factor Ce (say 1). 
b. Calculate the h, using equation (4.11). 


с. Calculate the tube wall temperature using the 
formula: 


d. Calculate a new correction factor, at the 
evaluated tube wall temperature 


e ОСЛА 
£ ` a ) 
w 


C 


63 











10. 


Tl. 


NOTE: 


E2. 


з, 


e. Repeat steps a through d until C, assumed 
in step a approximately equals Ce calcu- 
lated in step d. 


Calculate the log-mean-temperature difference, 
(LMTD). 


sat Тс. 


= - 12 
tn ( таак e) 2022) 


Calculate the overall heat-transfer coefficient 
using the formula: 


ы - а" 
О LMTD (2,13) 


Calculate the outside heat-transfer coefficient 
using the equation: 


p D Е (4 113) 
М 





The validity of equation (4.13) is based on the 
assumption of negligible water-side fouling 
resistance and the resistance due to noncondensables. 


Calculate the normalized, local, outside heat- 
transfer coefficient, ћу ћ) · 


Calculate the normalized average, local outside 
heat-transfer coefficient (бићу). 





С. DATA-REDUCTION PROGRAM 

A computer program was utilized to analyze the raw data. 
The program was in BASIC language and was run on an HP-9826 
computer system. A peripheral plotter was used to plot the ` 


results. The program is presented in Appendix D. 
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V. RESULTS AND DISCUSSION 


А.  SIEDER-TATE COEFFICIENTS FOR SMOOTH AND ROPED TUBES 

The modified Wilson Plots for smooth and roped tubes are 
shown in Fiqures 16 and 17. Of particular interest is the 
very Small intercept, which represents the estimated fouling 
factor. For example, the estimated fouling factor for data 
run STSD-11 (for smooth tubes) is 1.37x197© m¢-K/w 
based on the inside area. This value is only three percent 
of the typical value (4.4Х1075 т2.К/й) used for 
design purposes. This very small value supports the initial 
assumption of negligible fouling factor for the success of 
the modified Wilson-Plot method. 

Table III is a summary of the Sieder-Tate coefficients 
calculated for smooth and roped tubes. The average Sieder- 
Tate coefficient was calculated to be 9.929+9.081 for 
smooth tubes and 0.061-«0.002 for roped tubes. Thus. the 
roped tubes have a Sieder-Tate coefficient 2.1 times greater 
than that for the smooth tubes. This increase is mainly due 
to the increased surface area, turbulence and swirl effects. 

The Sieder-Tate coefficient derived for the smooth tubes 
is 7.4 percent greater than the value (0.027) published in 
the original. generalized correlation. This increase can be 


easily explained by the short condensing tubes used in this 
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TABLE III 
SUMMARY OF SIEDER-TATE COEFFICIENTS 
FOR SMOOTH AND ROPED TUBES 


| Sieder-Tate 
File Name Tube Type Coefficient 


STSD-2 Smooth 0.0287 
STSD-10 Smooth 070291 
STSD-11 Smooth 0.0296 
RTSD-3 Roped 0.0591 
RTSD-4 Roped 020589 
RTSD-6 Roped 0.0621 
RTSD-7 Roped 0.0627 
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experiment. Even though the flow condition was hydro- 
dynamically fully-developed, it was thermally developing 
throughout the condensing length, resulting in a greater 
Sieder-Tate coefficient. 

It is worth noting that the Sieder-Tate coefficient 
derived by using the original. Wilson-Plot method (i.e., 
plotting 1/0, versus l/X), consistently gave values about 
l0$ greater than the values obtained using the modified. 


Wilson-Plot method. 


B. SMOOTH TUBES 

Figure 18 shows the variation of the normalized, local, 
condensing coefficient for five tubes. The data points lie 
about 40 percent above the curve predicted by the Nusselt 
theory. This close agreement was considered to be an 
indication of the proper opepation of the test apparatus and 
the data reduction procedures. 

As discussed earlier in Chapter II. the Nusselt theory 
for a tube bundle is based on a number of basic assumptions. 
For example, in reality. the condensate does not fall as a 
continuous sheet as assumed for the Nusselt theory, but it 
forms drops before leaving the tube. This phenomenon has a 
tendency to result in a larger condensing coefficient than 


that predicted by the Nusselt theory. 
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Fiqure 19 shows the normalized, average condensing 
coefficient for Five tubes. йезе data points show a much 
smoother trend than that for local values. Further, these 
data points are closer to the Eissenberg correlation than 
the Nusselt prediction. 

Figure 20 shows the variation of the normalized, local, 
condensing coefficient for 30 tubes. The data points lie up 
to 35 percent above the curve predicted by the Nusselt 
theory. The outside heat-transfer coefficielt for the first 
tube was ten percent smaller than the value predicted by the 
Nusselt theory. 

Figure 21 displays the normalized, average condensing 
coefficient for 38 tubes. Again, these data points show a 
smoother trend than that for local values. These points lie 
between the Eissenberg and Nusselt predictions. 

Figure 22 shows the least-squares-fit curve for the data 
points for smooth tubes. The exponent derived is -0.154, 
and it is in close agreement with the value of -0.14 derived 


by Noftz [Ref. 12]. 


SEZZROPED TUBES 

Figure 23 displays the variation of the normalized, 
local, condensing coefficient for five tubes. The data 
points lie up to 15 percent above the curve representing the 


Nusselt theory. 
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Figure 24 shows the normalized, average, condensing 
coefficient for five tubes. These points show a trend very 
close to the Eissenberg relationship. 

Figure 25 represents the plot of the normalized, local, 
condensing coefficient under inundation conditions up to 30 
tubes. 

Figure 26 shows the normalized, average, condensing 
coefficient for 30 tubes. These data points lie about 
midway between the Eissenberg and Nusselt predictions. 

Figure 27 shows the least-squares-fit curve for the data 
points for roped tubes. Тһе exponent calculated is -0.183. 
Since the interest is in the large tube bundles. the curve 
fit was generated only for tubes 11 through 30. This is 
evident from the figure as the curve fit shows poor agree- 
ment when compared to the data points for the first ten 


tubes. 


D. ROPED TUBES WRAPPED WITH WIRE 

The outside. heat-transfer coefficient for the first 
tube was 11.3 percent greater than the value predicted by 
the Nusselt theory. This increase is in agreement with the 
manufacturer's claim that the special geometry of the roped 
tubes has a thinning effect on the condensate film, result- 


ing in a larger condensing coefficient. 
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Figure 28 shows the variation of the normalized, local, 
condensing coefficient for five tubes. The data points lie 
up to 63 percent above the curve predicted by the Nusselt 
theory. 

Figure 29 displays the normalized, average, condensing 
coefficient. The data points are well above the Eissenberg 
correlation. 

Figure 30 shows the normalized, local. condensing 
coefficient for a bundle of 30 tubes. The data points are 
Scattered within the limits of uncertainty, as predicted by 
the error analysis (see Appendix C). 

Figure 31 shows the variation of the data points 
representing the normalized, average, condensing coefficient 
under inundation up to 30 tubes. The data points lie up to 
100 percent above the curve predicted by the Nusselt theory. 

Figure 32 represents the least-squares-fit curve, which 
has a derived exponent of -0.039. The curve fit is in good 


agreement with all the data points. 


E. SMOOTH TUBES WRAPPED WITH WIRE 

The outside heat-transfer coefficient for the first tube 
was up to six percent greater than the value predicted by 
the Nusselt theory. This increase is caused by the thinning 
effect on the condensate film, resulting from the surface- 


tension forces acting toward the wrapped wire. 
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Figure 33 shows the variation of the normalized, local, 
condensing coefficient for five tubes. The data points lie 
up to 76 percent above the curve predicted by the Nusselt 
theory. 

Figure 34 shows the normalized, average condensing 
coefficient for five tubes. These data points are well 
above the Eissenberg correlation. 

Figure 35 shows the variation of the normalized, local, 
condensing coefficient for 30 tubes. The data points lie up 
to 107 percent above the curve predicted by the Nusselt 
theory. 

Figure 36 displays the normalized, average. condensing 
coefficient for 30 tubes. Тһе data points lie above the 
curve representing the Eissenberg correlation. 

Figure 37 shows that the curve fit is in good agreement 


with all the data points. The derived exponent is -0.037. 


F. OBSERVATIONS 
1. Smooth Tubes 

During all runs, complete film-wise condensation was 
observed without any visible evidence of drop-wise conden- 
sation. This was mainly achieved by the tube cleaning 
procedures used. 

A slow rate of condensate droplet migration. from 
cooling water outlet end to the inlet end, was observed at 


the bottom of each active tube. This problem was minimized 
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TABLE IV 


SUMMARY OF RUNS WITH INUNDATION 


External Wire 


File Name Dube. Type Wrapped or Not 


STNWI-1 Smooth No 
RTNWI-1 Roped No 
RTNWI-2 Roped No 
RTNWI-3 Roped No 
RTNWI-4 Roped No 
RTNWI-5 Roped No 
RTWI-1 Roped Yes 
RTWI-2 Roped Yes 
STWI-1 Smooth Yes 
STWI-2 Smooth Yes 
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TABLE V 


SUMMARY OF RUNS WITHOUT INUNDATION 


External Wire 


File Name Tube Type Wrapped or Not 
STNWEV-2 Smooth No 
STNWNI-1 Smooth No 
STNWNI-2 Smooth NO 
STNWNI-3 Smooth NO 
RTNWNI-3 Кореа NO 
RTWNI-1 Roped Yes 
RTWNI-2 Roped Yes 
RTWNI-3 Roped Yes 
STWNI-2 Smooth Yes 
STWNI-4 Smooth Yes 
STWNI-5 Smooth Yes 
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COMPARISON OF hy 


Tube Type 


Smooth 
Roped 
Roped 


Smooth 


TABLE X 


FOR UNINUNDATED TUBE RUNS 


Ди 


External Wire 
Wrap or not 


NO 
NO 
Yes 


Yes 
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О 156 
079019 
0790651 


0293359 








TABLE XI 


COMPARISON OF Ay FOR INUNDATION TUBE RUNS 
(Dy 
u 
Tube Type External Wire hy , 
Wrap or not hau 
Smooth No 059 
Roped NO 0.63 
Roped Yes 0.86 
Smooth Yes 0.86 


1:05 





TABLE XII 


EXPONENTS OF THE LEAST-SQUARES-FIT 


File Name Tube Type External Wire Exponent 
Wrap or not 

| STNWI-1 Smooth No 01153 
RTNWI-1 Roped No 0.183 
-2 Roped No 0.191 
-3 Кореа МО 0.179 
-4 Roped NO 0.195 
ктит-1 Roped Yes 0.039 
-2 Roped Yes 0.039 
STWI-1 Smooth Yes 0 309 7 
-2 Smooth Yes 0.056 
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after leveling the tubes, by adjusting the leveling nuts on 
the test condenser support bracket. When the condensation 
or the inundation rate increased, it was observed that the 
drops were formed at more sites along the tubes, but the 
droplet size for each tube was nearly the same. 
2. Roped Tubes 

Again» during the experiments, there was no evidence 
of drop-wise condensation. It was observed that the falling 
drops were coalesced on the lower surface of the tube, 
especially in the space between the two successive grooves. 
and then rivulets were formed and fell on the tube below. 
Under condensate inundation. the phenomenon of the rivulet 
formation was more intense. It is also worth noting that 
the droplet-formation frequency was higher near the cooling 
water inlet end. This can be easily explained by the 
larger, local, Sieder-Tate coefficient and by the larger 


temperature difference, T - фи at the inlet сопраг- 


sat 1 


ed to Т. - T at the outlet end. 


at со” 
3. Roped Tubes Wrapped with Wire 


The condensate was formed in the space between the 
Successive grooves wrapped with wire. Then due to the 
Surface tension forces. the condensate was drawn to the base 
of the wire, and drops were formed at the bottom surface of 


the tubes. 
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Under inundation conditions, the above-mentioned 
drop formation and movement were more intense. It is worth 
noting that no splashing was observed during the experiment. 

4. Smooth Tubes Wrapped with Wire 

It was observed that the condensate forms between 
two successive, helically-wrapped wires, and was then drawn 
immediately towards the space between the wire and the tube 
surface.  Rivulets were drawn from the base of the wire to 
the next tube below. It was observed that there was no drop 
migration along the tube. 

Under condensate inundation. the above mentioned 
drop formation was more intensive. but again no splashing 


was observed. 
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VI. CONCLUSIONS 


l. The average. outside. heat-transfer coefficient for 
30 smooth tubes was 59 percent of the Nusselt coefficient 
calculated for the first tube in the bank. 

2. The average, outside, heat-transfer coefficient for 
3@ smooth tubes wrapped with wire was 86 percent of the 
Nusselt coefficient calculated for the first tube in the 
bank. 

3. The average. outside. heat-transfer coefficient for 
38 roped tubes was 63 percent of the Nusselt coefficient 
calculated for the first tube in the bank. 

4. The average. outside, heat-transfer coefficient for 
30 roped tubes, wrapped with wire was 86 percent of the 
Nusselt coefficient calculated for the first tube in the 
bank. 

5. The Sieder-Tate coefficient for roped tubes was 2.1 
times greater than that for smooth tubes. 

6. Wire wrapping considerably improves the average 
condensing coefficient for both smooth and roped tubes in 
tube bundles. 

7. Of all cases investigated in this study, roped tubes 
with wire wrap would be the best candidate for deSigning 
compact condensers. However, the water side pumping power, 
which was not investigated in this work, could be consider- 


ably higher than that for the smooth-tube case. 
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VII. RECOMMENDATION 


A. TEST APPARATUS MODIFICATIONS 


The following test apparatus modifications are 


considered advisable: 


E. 


Redesign the test condenser flanges to allow a large 
sealing area, and also to allow for the possibility 
of further experiments with different S/D ratios. 


Redesign the existing test condenser hotwell to 
allow more reliable and convenient measurement of 
condensate flow rate. 


Install a larger heater for the perforated-tube 
water supply tank to facilitate the simulation of 
larger tube bundles. Also. modify the temperature 
control system for the perforated-tube water supply 
tank, to allow for more rapid heating and cooling to 
reduce long delays between runs. 


Install a glass window on the rear side of the test 
condenser, in order for the operator to have a view 
on both sides of the test tubes. 


Make the system vacuum tight so that data can be 
taken at vacuum conditions. 


B.  ADDITIONAL TESTS 


The following additional tests would be important in 


this continued investigation. 
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Conduct tests with enhanced tubes manufactured by 
Yorkshire Imperial Metals, Ltd. 


Take movies of the condensation process so that 
further conclusions could be drawn with regard to 
condensate drop phenomena and their relationship to 
condenser performance. 


Conduct tests varying the pitch of the wrapped wire 


on smooth titanium tubes and determine the optimum 
wire pitch. 
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Withdraw gas samples from the test condenser and 
analyze on the gas chromatograph to determine the 
effect of noncondensable gases on the condenser 
performance. 


Conduct tests to investigate the effect of vapor 
velocity on the heat-transfer coefficient. 
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APPENDIX A 


OPERATING PROCEDURES 


INITIAL PROCEDURES 


I. 


Energize the main circuit breaker located in Power 
Panel P-2 on the wall to the right of the test 
apparatus. 


Energize the circuit breaker on the left side of the 
old control board by pressing the ON button. 


Energize the following switches in the control 
panel: 


а. #1 - Perforated tube water supply tank (feed 
pump). 


b. $2 - Outlets. 

с. #3 - Hot water heater. 

d. #4 - Condensate pump. 

e. #6 - Cooling tower. 

f. #7 - Cooling water pump. 

Ensure all test apparatus valves are closed. 

Fill the perforated tube condensate supply tank with 
distilled water. Set the temperature controller for 
the perforated tube water supply tank at about 

959 С, fully open the recirculation valve, P-l, 

and start the pump to begin heating the water. The 
controller will have to be reset to the proper 
supply temperature once steady-state conditions are 
obtained. 


Fill the cooling water supply tank. This can be 
done by backfilling with valves CK-1 and CW-4 open. 


Energize the data acquisition system. 


daly 





B. OPERATION 


M 





House Steam 
a. Open the main supply valve. 


b. Open valve MS-3 until the pressure gage 
indicates the desired steam supply pressure. 


с. Fully open MSD-1 to drain any condensate. 


d. Open valves MS-4 and MS-5 until the desired 
Steam supply pressure is obtained and re-adjust 
MS-3 as necessary. 


Condensate System. To collect the condensate in the 
test condenser hotwell, operate the system with 
valve C-1 closed. After a test run is completed, 
open valve C-1 to drain the condensate by opening 
valve C-2 to the bilges. or operating the condensate 
pump in order to fill the perforated tube supply 
tank. 


Cooling-Water System 

a. Open valves CW-1. CW-2 and CW-3. 

b. Ensure valves СК1-1 and CW-4 are closed. 
c. Energize the two cooling water pumps. 


d. Open valves CW-5, CW-6. CW-7. CW-8 and CW-9 to 
obtain the desired cooling rates. 


Perforated Tube Water Supply System 


a. Once steady-state conditions are achieved, reset 
the temperature controller to the proper 
inundation temperature. 


b. Adjust the rotameter to the required flow rate 
for each run. The supply tank recirculation 
valve may have to be adjusted to achieve the 
desired flow rate, but should never be fully 
closed to avoid damaging the pump and to ensure 
uniform water temperature in the tank. 


c. Refill the supply tank as required. by using the 
existing piping filling system, by operating the 
condensate pump; or filling the supply tank with 
distilled water filling line. 
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Miscellaneous. 


To maintain a clear test condenser window, open 
valve A-2 and then energize and adjust the air 
heater power supply. When securing, always turn off 
the power supply first and allow the air heater to 
cool before securing valve A-2. 


SECURING THE TEST APPARATUS 


Ш. 


10. 


ІНІ 


Secure the steam valves MS-5. MS-4. MS-3 and the 
main supply valve. 


Secure the air compressor. 

Secure the perforated-tube water supply system by 
securing the pump, temperature controller, and 
valves P-l and P-4. Drain the supply tank by open- 
ing valve P-3 (if desired). 


Secure the test condenser viewing window air heater 
as prescribed above. 


Secure the data acquisition system. 

Allow the test condenser to cool down for about 30 
minutes. then secure the cooling water pumps and 
close valves CS-1. CW-2, CW-3. CW-5. CW-6, CW-7,. 
CW-8, CW-9 and C-4. 

Drain the test condenser hotwell. 

Secure all circuit breakers. 

Drain the cooling water system piping and rotameters 
by opening valve CW-10 and leaving open the cooling 
water rotameter supply valves. 


Drain the cooling water supply tank by opening the 
drain valve via the remote operating rod. 


Ensure all valves are secured. 
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APPENDIX B 


SAMPLE CALCULATIONS 


A. RUN STNWNI-1 


A sample calculation is performed in this section to 


illustrate the solution procedure used in the data reduction 


program [Ref 29]. 


The STNWNI-1 run was selected to perform this analysis: 


INPUT PARAMETERS 


File 

Pressure condition 

Inundation condition 

Month, date and time 

Run number 

Tube number l 2 
Inlet temp (9C) 28.53 28.54 
Outlet temp (°C) 32.55 32.20 
Saturation temperature 

Degree of superheat 

Condensate temperature 

Static pressure 


STNWNI-1 
Atmospheric 
5 tubes 
05:11:10:08:50 
l 
3 4 5 
28—55 28.54 28:58 
92.17 32205 319583 
100.24 (°c) 
1.88 (ОС) 
92.41 (ОС) 
767.33 (mm Hg) 


The following calculations are further limited only to 


the first tube. 


ресе пајтоп о Average Bulk Temperature 


T, (1) = Тс; (1) + Тео (1) 
T, (1) = (28,53 + 32.55) 0.5 
T, (1) = 30.549C 


115 


x0.5 





2. “Шегпеоризезса1і Properties 


Р = 5.329 
r 
3 
p = 995.2 kg/m 
Io oS M SUA 
E = 4.178 Kj/kg.K 
3 


k, 7 619 x 10 7 W/m-K 


NOTE: All properties are calculated at the average bulk water 
temperature from Table A.6 p. 782 [Ref 30] 


3. Cooling Water Mass Flow Rate 


m = 14.22 Kg/min = 0.237 Kg/s 


4. Determination of Cooling Water Velocity 


V ще 
w PA. 
1 
1 
(14.22) 60 
М -4 
(995.2) (1:56 x 10 `) 
үле = 1.53 m/s 


5. Determination of Reynolds Number 


aD: 
R ee = 
e m 
8 = (995.2) (1.53) (0.141) 
e 798 x 10`? 


116 





R_ = 26,904 
E 


Determination of Heat Transfer 


SM 5 Ті) "Соч 
О = (14.22) Е (55 28.53) (4.179) 
Q 2 3,980.5 W 


Determination of Heat Flux 


"= 
а т ет, 
О 
· – _____3,980.5 ____ 
а ms 0.015875) (0. 05) 
q" = 261,682.2 — 
m 


Determination of Nusselt Coefficient 


Dd Dorem +з 
å Е £ £g 
Buu TES 0.651 u D e 
£ O° 
Assume Te = Ds 
т; = 100.249с 


р; = 957,8 kg/m" 


_ 3 
he, = 2.2564 x 10” J/kg 


Ke = 682 x 107? W/m-K 


1257 





А 1076 М. 
2 
m 


NOTE: All properties are calculated at the Saturation 
Temperature. From Table A.6 p. 782 [Ref 30] 


3 2 123 


3) - (957.8) (2256.4х103) (9.81) 


Be = 0.651 (682x10 ч 
(281x10 9) (0.015975) (261682.2 
hy,  11492.6 W/m?k 


0, Determination of Te E 


— = E 

Tg c ^ Tsat hu шы 
a _ 261685.3 
= O 

ТЕ e ӨЗ. 89°C 


10. Thermophysical Properties 


ke= 673 = 1077 мутк 


О = 963 Kg/m? 

u = 326 x 10 9 N.s/n* 

h. = 2.289.5 x 10° J/kg 
E .289. 


NOTE: All properties are calculated at the film temperature. 
From Table A.6 p. 782 [Ref 30] 
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B. 


ие. 


№». 


ша. 


Determination of Nusselt Coefficient 


3.2 1/3 
Не св 5 *5929.081 
Nu He-Dg-q 
MU (ог вз) (2259:5-х 107) (981) 
ив (326 x 10 9) (0.015875) (261685.4) 
2 

hj, = 10,885.3 W/m^.K 
Determination of ТЕ с 

E: pw 
Те с B Ec h 0.5 

Nu 

B _ 261685.4 
ТЕ с = 100.24 l0888.3 9-2 
T. , - 88.229C 

Pc 


Determination of Logarithmic Meal Temperature 
Difference (LMTD) 


Teo Т Tei 
ш! sat сі) 
Teat ‘co 
_ 32.55 = 28.53 
ыы zu (100.24 - 28.53, 
100.24 = 32.55 
LMTD - 69.689C 


Determination of Overall Heat - Transfer Coefficient 


Чо = MTD 
y = 261685.4 
о 69.68 


о 





5. 


I. 


I. 


19, 


y, = 3755.5 —— 
m K 


Determination of Inside Heat-Transfer Coefficient 


Assume Ce = 1.1 


С; = 0.029 
m О 
h; = 5. Са Ке P. "СЕ 
T deu 0.8 0935399. 
h, = 92012 (0.029) (26904 °°°) (5.829 mm 
h, = 9303.2 W/m^K 


Determination of Inner Wall Temperature 


c 

II 

H 

+ 
үа. 

| 


(261685.4) (0.015875) 
w ^ 30.54 + —79303.2) (0.0141) 


= 
|| 


Т 52:229C 


w 


Determination of u,, at the Average Wall Temperature 


6 


um (62.20C) = 453 x 10 ТТТ 


Determination of Correction Factor 


и 0.14 
С = W : 
Бе 5 К 
: w 


798 x 10 9] 9:14 
Ecc -6 

453 x 10 
DE EIOS 
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19. Determination of Outside Heat-Transfer Coefficient 


1 
МЕГЕ 





- "EE Е 
UG D. E w 
Ш = i WEE 
3785-8 - T0-1813313393-27 - 0-000042925 
h = 9772.3 W/a^k 


В. RUN STSD-11 


A sample calculation is performed in this section to 


plot program [Ref 29]. 
The DP-11 run was selected to perform this analysis 


INPUT PARAMETERS 


File DP-11 

Month, date and time 04:01:13:57:20 
Data point #1 

Steam Saturation Temperature i 100701 C 

Inlet Temperature 24.420C 

Outlet Temperature 30.759C 
Flowmeter Reading 10$ 


1. Determiantion of Average Bulk Water Temperature 


T = (T . +T ) 0.5 
C 


b cu О 


(2442. E 30. O o 5.0 27. 5 9OC 


In 


2. -Themmophysical Properties 


Ва - 5.83 


K 
р + 997 = 
m 


illustrate the solution procedure used in the modified Wilson 
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Coy = 4.180 Kj/Kg-K 


= E E. s/n 


k = 0.613 W/m-K 
М 


NOTE: All properties are calculated at the average bulk water 
temperature from Table A.6 p. 782 [Ref. 30] 


3. Cooling Water Mass Flow Rate 


М. = 6.36 kg/min 


£ 
4. Determination of Cooling Water Velocity 
и на 
Мм PA, 
5 
1 
(6.36) 60 
КҮТ mI ———41 
М (997) (1.56 х 10 
s. - 0.68 m/s 


5. Determination of Reynolds Number 


Рио: 
Ке E но 
W 
a = (997) (0.68) (0.0141) 
е SE TE 
A 
e 


6. Determination of Heat Transfer 


Q = Me у Bec y Lu ` Cow 

ие №. 36) (=) (30.75 - 24.42) (4180) 
Q 

О = 2,805 И 


4522 








t9. 


Determination of Logarithmic Mean Temperature 
Рекао БИО 


ША =.=-Т 


LMTD - — 
= sat сі) 
Tsat Тсо 
Же 815-242 


n Ga 100.01 - 24.42 
POCO = 30.75 


LMTD = 72.389C 


Determination of Overall Heat - Transfer Coefficient 


U ee ae a. 
О T-DG: b: LMID 


_ 2805 
о” т(0.015875) (0.305) (72.38) 


E W 
UG Z9 0477415 о 


Determination of Sieder - Tate Parameter 


хор -0.8.р -0.3333 
e E 

 - (11/154:0-9) (3.8377 >) 

x = 0.0003212 


Determination of Inside Heat - Transfer Coefficient, 
based on assumed Sieder - Tate Coefficient 


Assume Ce = 1.1 апа 
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EL. 


ша. 


ES. 


14. 


-3 
О ОЗ 1540-38) (5.330-333 


; 7 4,463.4 > 
m“ ° K 


a 
| 





Determination of Average Wall Temperature 


E Q 

Ta Е Тр с пе. 
1 l 
a 2805 
М b п0.0141) (0.305) (4,453.4) 
= = О 

Т. T, 46.519C 
T = 46.51 + Т, 
T = 46.51 + 27.59 = 74°C 


Obtain uu at the Average Inner Wall Temperature 


к он 2 
UE CD) 2 37572x- 10 N-s/m 


Determination of Cec 


uy \ 0.14 
В“ 
Ес ( (T > 
W Ww 


А (г: x ш) 0.12 
E Е -6 
37.522510 


19102267 


Сес 


Iterate for h; 


Assume Ce - 1.12 апа 


er 0.03 


T: 
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Ше. 


6. 


7. 


Ku 0.8, 0.3333 
=" с, | 


1 . 1 e t f 
1 
es 613 x 10° (0.03) (11,1549-®) (5.830-2333) (1.12) 
4 0.0141 . r . : 
КИШ Бе 97 w/m K 


Determination of Wall Thermal Resistance, Based on 
the Outside Diameter 


(D -D.) D 
Da . o O 
wW к .(р -р.) 

m О 1 


(0.015875 - 0.0141)(0.015875) 





ЦИ ra oo. 015875 + 0.0141) - 
m^'K 

В = 0.000042925 

w W 

Assume R = O 


£, 


Determination of Outside Heat - Transfer Coefficient 








h = 1 
D 
і C 2 
U w D. bs 
Q PP] 
пије 1 
нысы = bU Оп 
acm Scc OT (4163.4) 
ПИ 273 вик И 
5 
m -K 


Set Qo = Q = 2805 W 


Determination of Actual Sieder-Tate Parameter 


X 


Х = ---- 


Cte 


00009212 
1.12267 
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x = 0.0002861 


INPUT PARAMETERS 


File STSD-11 
Data point #2 
Steam saturation temperature 100.03°С 
Inlet temperature 24.390C 
Outlet temperature 29.460C 
Rotameter setting 15$ 
18. Determination of Average Bulk Water Temperature 
= . + Š 
T. (т, Тсо) 2085 
Ты - (24,39 ғ 29,46) 0.5 
= O 
T, 26.920C 
19. Thermophysical Properties 
P_ = 5,823 
r 
E 3 
и = 996 Kg/m 
Cow = 4.180 Kj/KgK 
Der ier” N.S/m^ 
к, = 613 х 1073 нупк 
20. Cooling Water Mass Flow Rate 
Me = 9.72 Kg/min 
22. Determination of Cooling Water Velocity 
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22. 


28. 


24. 


25. 


I 
(9.72) 60 


Vo Т -4 
(996) (1.56 x 10 ^) 
Мы = 1.04 m/s 


Determination of Ro 


P V D. 


R. = = i 
w 

R = (996) (1.04) (0.0141) 

= 855 x 10 9 

R, = 17,082 


Determination of Heat - Transfer 


Q = Me у ES Те; ) Cow 
Q = (9.72) (5) (29.46 - 23.39) (4180) 
Q 2 3,433.2 W 


Determination of Logarithmic Mean Difference 


Temperature (LMTD) 


Teo ^ tci 
LMTD = — 
га ( sat “cl ) 
T -T 
sat “со 
š 29:46 - 21-6 
кр = үп(100.03 24.39 
100.03 - 29.46 
LMTD = 73.08°C 


Determination of Overall Heat-Transfer Coefficient 


A 
O T- DoE. HMTD 
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26. 


27. 


28. 


29. 


0995-2 
СОТО БТ) (00870757) (73.08) 


U 3,088.4 W/m^K 


O 


Determination of Sieder-Tate Parameter 


08 20895353 


Х = Ro P. 
x = (17,082 0:9) (5.823 03333, 
X = 0.000228 
Determination of = 

i 

1/3 

т. ES - (25) 
h; UG h. Q w Do 
[Dm >. | ЩЕ " Т ( 2805 2400004295 | 0.015875 
h; Б 3088.4 10273 \37 133.2 i 0.0141 
1 m*K 


Determination of Average Wall Temperature 


Q 


п = 
w b 7.D,-L-h, 

3433.2 
ВЕ 20:92 09015975): (00:305) 01,77 322€) 
T = 74.200C 
w 


Determination of и „а the Average Wall Temperature 


-6 2 
ACT) = 375,2 x 10 N-s/m 
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50 


За. 


Determination of Cfo 


u 0.14 
Ф = И. СУВЕ 
s Ur Tw) 
NEG 0089 0.14 
С с шш: 
375.2 x 10 
rs 07 


Determination of Actual Sieder-Tate Parameter 


x = g— 
ЕС 
х - 0-000228 
1 124159 
X - 0.0002027 


T0209 


E T ||" 


P, 








APPENDIX C 


UNCERTAINTY ANALYSIS 


The general form of Kline and McClintock [Ref. 31] "second 
order" equation is used to compute the probable uncertainty in 


the results. For some resultant, R, which is a function of 


primary variables Xi, Хор... X the probable uncertainty in 
КОШСОК is given by: 
122 
ек 22 eR 2 OR 2 
= | (= a о 
SR | 9X, ) + ex, 8X] +... + (5х Xa? | (Са) 
n 

where 25 5Х.,..., 2 is the possible uncertainty in each of 


the measured variables. 


A. UNCERTAINTY IN THE COOLING WATER VELOCITY 


| m 
E PA. 
ji 


Applying equation (C.1) the following equation results: 


2 


ô V ° o A: 
м _ ôm, 2 TN jiu 
КА = ГЕ E (—с) F E | (С.2) 


The following uncertainties were assigned to the variables: 


ám = 40.01 kg/s 
=> kg/m? 
SA, = +0.0001 mf 
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B. UNCERTAINTY IN THE REYNOLD'S NUMBER 


_ 1 
e Hu 


R 


The probable uncertainty is given by: 


172 
SR SV Ó DE 
5 
ER = Qe + [e + (=r = Put (е 
р V D. 
e w a 
The following uncertainties were assigned to the variables: 
ee +3 kg/m” 
50; = OOOO En 
бы = +8х107® N-S/m* 
Bs = from equation (C.2) 
C. UNCERTAINTY IN HEAT TRANSFER 
Ола Тео Тај) Сри 
The probable uncertainty is given by: 
p» 
: óT ПЛ SC 
i0 6m, 2 со 2 Ст? м, 2 
us [ E) + (тзт) + (т ст -) EE C ) | (C.4) 
CONS Go ci pw 


The following uncertainties were assigned to the variables: 


ám = +0.01 kg/s 


= о 
Se 9025 “с 
ЕЕ! = 025 °C 

cu E 
$C 7 t8 J/Kg'C 
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The 


The 


The 


The 


UNSBRTATNEY OF THE HEAT-FLUX 





probable uncertainty is given by: 
lao 
за" 211 (40,2, (“о)2, (8,2 (C.5) 
e m Q D. L 
following uncertainties were assigned to the variables: 


6D. = +0.0001 m 
Q == 


61 


+0.0001 m 


6Q as found from equation (C.4) 


UNCERTAINTY OF h 
Nu 


ЗА "а TAS 
h 2 0.65] | ———ÓÉ— 
Nu Н 0 9 


probable uncertainty is given by: 








h k 6h Su 
NU = 0.651 сто)? + (2 592 + (L-£2)* . (88)? + (FE)? 
Nu f fg Hr 

9 Do 2 ба" 2 - 


following uncertainties were assigned to the variables: 
6ke = +0.0012 w/m’k 
бр = +3 kg/m? 


ô h 


+0.48 J/kg 
6 


fg 


Su, = +8x10°° w-S/m' 


6D, = +0.0001 m 
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The 


The 


The 


The 


5g = 40.001 m/s? 


óq as found from equation (C.5) 


UNCERTAINTY OR T.. 
film, 


f sat 


ая 
T = T - — 0.5 
С huu 


probable uncertainty is given by: 


1/2 
Е ST Е sh 
= = Ë + (292) + un? | (cam 
в sat d Nu 


ST 





T 
C 


following uncertainties were assigned to the variables: 


ST = +0.025 °C 
sat -- 
6q"- as found from equation (C.5) 
Shu = as found from equation (C.6) 


UNCERTAINTY OF OVERALL HEAT-TRANSFER COEFFICIENT 


Чо = LMTD 


probable uncertainty is given by: 
1/2 
806 6q",2, ,§(LMTD) ,2 
— = ( а") + (>=) (C.B) 


Чо LMTD 


following uncertainties were assigned to the variables: 


óq = as found from equation (C.5) 


$ ( LMTD) as found from equation (C.9) 
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H. UNCERTAINTY FOR LOGARITHMIC MEAN TEMPERATURE DIFFERENCE 


( LMTD) 
Bone ci 
LMTD = тт 
ар ) 
S “СО 


The probable uncertainty is given by: 


SLMTD _ ( uere x) ) E 
LMTD ПО Со Ab EU 
In ! 
S "co 
2 
Ты 
i (m em T -T ) B 
s ci у (5 сі, 
Mp =p 
6 “СО 
1/2 
уе ^ 
E e -T _) T -T . ) (С.9) 
в co (S ci, 
In? m 
Saco 


The following uncertainties were assigned to the variables: 


6T. — +0.025 “С 
6Т_. = +0.025 °С 
C 1 T 
ST = +0.025 °C 
CO Би 


ШІ ГГЧСЕКТАТМТҮ OF INSIDE HEAT-TRANSFER COEFFICIENT 


Dow CR 


T De T ч 


kw 0.8 
e D 


0.333, 5,0.14 
uw 
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The probable uncertainty is given by: 


The 


The 


gh. Š k ôD. 0.858 0.3336Р $C. 
ү м, 2 Ш е, 2 Ба 2 
À -— E DENS ie aj ae + (—— 5) «+ (в) + 
T ү J. e r т 
1/2 
0.146(u/u_) 
ыы 2 (C.10) 
W/E 
following uncertainties were applied to the variables: 
єк, = +0.0012 w/m*k 
SD; = +0.0001 m 
6С; = +0 0001 
SR, = as found from equation (C.3) 
or = ml 7 
B EID 2 
6 (м/н, ,) = 8х10 Ns/m 
UNCERTAINTY IN TEMPERATURE DIFFERENCE 
qe D. 
E n D. 
т I 
192 
Sh. ӛр ӛр. 
SAT _ да", 2 12 GI que 
AT = [ а" 15 n З E De | (ИЕ 


following uncertainties were assigned to the variables: 


öq"= as found from equation (C.5) 


sh, = as found from equation (C.10) 
80 = +0.0001 m 

O — 
60. = +0.0001 м 


I 
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The 


UNCERTAINTY IN OUTSIDE HEAT-TRANSFER COEFFICIENT 














О 1 D 
m = OR к 
EC D.h ы 
Iel 
SAL 59. 2 
NEC ( 2 ) P 
О UL (i _ R - ) 
U ү D.h 
О Tung 
SR. 
. ( 572 
Ш.с 
U ү D.h 
О ШП 
1/2 
D. oh 2 
ary M 
T 5 ic. 129 
ee ии E 
U Мм р.п. 
О J L 


following uncertainties were assigned to the variables: 


50 = as found from equation (C.8) 
sh, = as found from equation (C.10) 
SR, = +0.00001 m“ -k/w 


UNCERTAINTIES FOR THE NORMALIZED LOCAL HEAT-TRANSFER 


COPPPTICIENT h/h; 


This ratio is simply the heat transfer coefficient of a 


given tube, N, divided by that of the first tube, i.e., for 


the fifth tube, N=5 and: 


Ay 55 


h h 
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An application of equation (C.1) results in the following 


equation: 
6 (hy /h; ) sh, sh nae 
N^ 1 12 112 
hy/ Pi i Ay 
Note: Equation (C.13) is valid only for N7/2. For example: 
112 
Anh m | E. 24] 
һ,/һу h; h, 


M. UNCERTAINTY IN THE NORMALIZED AVERAGE HEAT-TRANSEER 

COEFFICIENT, h,/h, 

The normalized average heat-transfer coefficient is obtained 
for the Nth tube by taking the average of the heat-transfer 
coefficients of the first N tubes and dividing this by the 
heat-transfer coefficient of the first tube: 


Ay _ (hy + h, T... +t hy) N 


h h 


ЕН 


р 


Applying equation (C.1) to the above, the following equation 








results: 
Е N = 1/2 
б/р) Shi Sh. 2 
т B 2 -N (Ен) (C.14) 
ma)! ігі yh $ 
i 1. 
1-1 
where N - the tube number. For example: 
B 72 
uA _ | Bi ou. q? | 
(Hy /hy) a hy +h, 21 
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For Run STNWNI-1 


V, = 1.54 + 0.11 "/s 


Re 


27546 % 2121 
О = 3934 + 169 м 


q" - 267830 + 11516 "/ 2 


Е W 
hau = 10802 + 164 / mek 


LMTD = 69.7 + 0.6 °C 


с 
|| 


3843.7 + 165 


W 
DEEP ова / 2. 


c 
|| 


от 
|| 


W 
11361.4 * 1249 "/ 2, 


~ 
|| 


0.8891 + 0.14 


~ 
|| 


0.9446 


|+ 


0512 
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APPENDIX D 


COMPUTER PROGRAMS 


A. DATA REDUCTION PROGRAM 


1000! FILE NAME: DRP 

1010! REVISED: May 20. 1983 

IU CDM /Ci/ CC7) 

1030 DIM тс: (2) „Тсо(4,2),171 4), МР+( 4), Ми( 4), НоС4) 

а рМ То са), nc». Тедо R3 (4), R4(4) ,S2(4) ,S4(4) 

1050! 

1060! ASSIGN COEFFICIENTS FOR THE 3-TH ORDER 

1070! POLYNOMIAL FOR TYPE-T «COPPER-CONSTANTAN) 

1080! THERMOCOUPLES 

1090 DATA 0.10086091,25727.94369.-767345.3295.78025595.31 

1100 DATA -9247485589,6.97688Е+11,-2.66192Е+13,3. 940722 +14 
1110 READ C(+) 

1120! 

ПАЗ ТОМ РОСЕ-оСАКЕ FLOw RATES TuROUCH THE 5 

1140! FLOH METERS (ka/min) 

[150 DATA 66.26,73.25.72.44.72.52,72.24 

1160 READ Mft(») 

11270! 

1180! ASSIGN SIEDER-TATE COEFFICIENT AND EXPONENT: 

1190! FOR REYNGLDS NUMBER 

1200 Ci=.029 

1210 Ex=.8 

1220! 

1230! ASSIGN GEOMETRIC VARIABLES 

1240 Di=.0141 ' Inner diameter (m) 

1250 Do=.0!5875 Outer diameter (a) 

1260 Ktn=21.3 ' Thermal conductivity of titanium (W/m-K) 
1270 | =.305 ' Condensing length (m) 

1250 Мс-3.3333 ! Number of unit cells across condenser width 
1500 ve 1. 5 ‘ Transverse tube pitch-to-diameter ratio 

1 ! 

1310! COMPUTE THE MINIMUM STEAM FLOW AREA IN THE TEST CONDENSER «м `2) 
1320 Amf=Nc*=Do*=(Pt-PI/(4*Pt)>>=L 

13908 

1340! COMPUTE INSIDE AREA AND WALI. RESISTANCE 

1350 Ai=PI*D1*2/4 

1360 Rw=DoxLUG(Do/Di)/(2*Ktm) 

1370! 

1386 PRINTER IS 703 

1390 CLEAR 709 

1400 BEEP 

1410 INPUT “ENTER MONTH, DATE, AND TIME (MM:DD:HH:MM:SS)",TimeS 
фа UEEDUIPUL 709:"TD'";TimeS 

1430 BEE? 

1440 INPUT "ENTER THE INPUT МОРЕ (1=3054AA,2=FILEI". Im 

1450 IF Im=2 THEN 

1460 BEEP 

1470 INPUT “ENTER THE NAME OF THE EXISTING DATA FILE" .OlddataS 
1480 PRINT USING ~10X.° "This analysis was performed for data stored in file "", 
104” :Dlddatas 

1490 ASSIGN AF¡le2 TO Olddatas 

1500 END IF 

1510 IF Im=1 THEN 

1520  BEEP 

1530 INPUT "GIVE A NAME FOR THE DATA FILE TO BE CREATED", MewdataS 
1540 CREATE BDAT tlewdatas.20 

1550 ASSIGN 3Filel TO Newaatas 


ша 





1560 
1570 
1530 


1680 
1690 


1720 
1730 
1740 
1750 
1760 
1770 


1800 


3U TUBES 
3) TUBES) 


END IF 

BEEP 

en Gives ne РОК THE DUTPUT FILE”, File_outS 

INPUT "ENTER THE PRESSURE CONDITION ¢!=47TM.2=VACUUM)” . Мр 
IF Mo*s1 THEN PRINT " Pressure condition: ATMOSPHERIC” 
= ПЕВ THEN PRINT " Pressure condition: VACUUM" 

EP 

INPUT “ENTER THE INUNDATION CONDITION (1*5 TUBES. 2-30 TUBESD".Mi 
IF M12 THEN PRINT " Inundation condition: 

IF Mi*1 THEN PRINT " Inundation condition: 
CREATE BDAT Fiie_outS.6 

ASSIGN »File3 TO File_outS 

Ja=0 

Мгип=0 

FOR I=0 TO 4 

S$3¢1)=0. 

S4¢i =u 

NEXT I 
Repeat: ! 

Mrun=Nrun+1 

OUTPUT 709: TD" 

ENTER 709; Times 

PRINT '" " 

PRINT USING "10X.,'""Month, date. and time: "",1580":TimeS 
IF Im=2 THEN Rdf 

BEEP 

INPUT “ENTER FLOW METER READINGS (AS PERCENTAGES)” .Fnt.Fm2 
run MOD 5=1 AND Mi=2 AND Nrun?5 THEN 

EOUM ENTER FLOW RATE FOR POROUS TUBE (AS A PERCENTI",Fpt 
QUTPUT wFilel:Fpt 

Mpt=-3.361613+10.076742*Fpt 

END IF 

DISP “START COLLECTING CONDENSATE” 

BEEP 

WAIT 20 

ПОТРИТ 709:"AR AFO ALI9" 

ШОР ЕТІ RI T1 Z1 FL1" 

READ INLET WATER TEMPERATURES 

FOR I=0 TO 2 

OUTPUT 709:"AS SA” 

EMER 722:Тс1 (1) 

CALL Pvsv(TeicT)) 

Тс: (Т) -МТепр <Тс: (Т),Т) 

NEXT I | 

READ OUTLET WATER TEMPERATURES 

Iı=2 

FOR I=0 TO 4 

IF I=0 OR I=32 THEN 

lu=2 

BLSE 

Iu=1 

END IF 

FOR J=0 TO lu 

Iı-Ilit! 

ООТРЏТ 709; "AS SA" 
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2170 ENTER 722:Tco(I.J) 

2180 CALL Tvsv(CTcoCI,J)) 

2190 Tco(I,J»?-FNTemp(TcoCI, AD, Ii) 
2200 NEXT J 

2210 NEXT I 


2230! READ STEAM TEMPERATURES 


МОШЕ РОК 1=15 ТО 16 

2260 ООТРОТ 709;"AS SA" 

Bae CENTER 722:Ts€(1-15) 

2280 CALL Tvsv(TstI-15)) 
2230 Ts(tI-15)=FNTemp(TstI-15),D 
2300 NEXT I 


320! READ CONDENSATE TEMPERATURE 


2340 OUTPUT 709:"AS SA" 
2350 ENTER 722:Tcon 

2360 CALL Tvsv(Tcon) 

2370 Tcon=FhNTemp(Tcon, 1?) 


2390! READ VAPOR TEMPERATURE 


2410 OUTPUT 709:"AS SA” 
2420 ENTER 722: Tv 

2430 CALL Tusv(Tv) 

2440 Tv=FNTemp(Tv,18) 


2460! READ VAPOR PRESSURE 


2480 DUTPUT 709:"AS SA" 
2490 ENTER 722:P_volts 


2510! COMPUTE AVERAGE WATER TEMPERATURES AT INLET 


2530 Tı(0)=-Tci(0) 
2540 ІТІС12-(Тс1(09)%Тс1(1929%.5 
ЙӘӘШ Т1(29-Іс1(1) 
eee 1i03)=(Feiti)+Tci(2))*.5 
2570 Tı(4)=Tcıt2) 


2530! COMPUTE AVERAGE WATER TEMPERATURES RT OUTLET 


2610 FOR 1=9 TO 4 

2620 IF І-0 OR I=3 THEN 

БН ток Тсо(Т.0) +Тсо(Т, 1) +Тсо(Т,2))«. 3333 
2640 ELSE 

2650 To(I)=(Tco(I,O)+Tco(1I,1))*.S 
2680 END IF 

2670 NEXT I 

2680 Tsa=(Ts(0)+T=(1))=.5 

2690 Pvap-FNPvsv(P volts) 

2700 Tsat=FNIvsp(Pvap) 

2710 Dsup-Tv-Tsat 


2730! READ INFORMATION FOR CONDENSATE FLOW RATE 


ӘТЕШ BEEP 
2 ogee EU t ENTER INITIAL AND FINAL LEVELS IN HOT WELL 1'".H1.H2 
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2770 Dh=H2-H1 

2780 IF Nrun MOD 5*1 THEN Msum*Ü 

2790 Hf1s540.4826*Dh 

2800 Mdl=Mf!-FNRhow(Tsat-10)#1.0E-6/60 

2810 Msun=-Msum+Mf1 

2820 IF Mi=2 AND Nrun<>30 AND Nrun MOD S=0 THEN 
2830 Mave=Msum/S 

2840 Set=(Mave*FNRhow(Tsat-10)/10"6+.03238)/.042132 
2850 END IF 


2860! 

2870 Rdf: ! 

2880! 

2890 PRINT USIHG "10X,""Run number = "'",DD':Nrun 

2900 PRINT " Tube + : 1 2 3 4 Se 


2910 If Im=2 THEN 

2920 IF Nrun MOD S=* AND Mi=2 AND Nrun>5 THEN ENTER ЗЕ 1 | ед; го 

Ж ЕЕЕ ЕНТЕН @Fiiec:1i1<(*).Tot=).isa.icon, iv. fvap,isat,Dsum,Fmi .Fa2 

2520 ENTER 3F1le2:H1,H2 

2950 END IF 

2950 IF Iim=! THEN QUTPUT ƏFilel;Ti(w),To(*w*),Tsa,Tcon,Tv.Pvap.Tsat.Dsup.Fml,Fm2 
КОО PRINT USING "70%, ""Inlet temp (Deg C?) :””,S5(DDD.DD.2X>":Tit=) 

feo PRINT USING "10X." "Outlet temp (Deg C):"",.5¢(¢0DD.D0D.2K)";: Tot =) 

2990 PRINT USING "10X.,""Saturation temperature cnp "t" Срва СУ Теа 
3000 PRINT USING "10X,""Degree of superheat “",20.00."“ (Dea C)""":Dsup 
3010 PRINT USING "10X."" Condensate temperature МЛ” (Пес Се Тер 
3020 PRINT USING "10X,""Statıc pressure ~" 3D DD. Snn Hg) ЕРуар 
3030 IF Im=1 THEN OUTPUT &F11e1:H1.H2 


3050! CALCULATE AVERAGE BULK TEMPERATURES 


3070 FOR I=0 TO 4 

3080 Tb(I)=(Ti(I)+To(I))*.S5 

3090 NEXT I 

3100! 

3110 IF Mi=! OR (Mi=2 AND Nrun<6) THEN As!=0. 
3120 IF Mı=2 AND Nrun>S THEN S1=As1 

3130 S!=As! 

3140 FOR J=0 TO 4 

3150 IF J=0 THEN Cuf=Fm1 

3160 IF J=1 THEN Cuf=Fme 

3170 Hf=Mft(J)=Couof /(100=60) 

3180 Tx=Tb(J) 

3190 Yiuw(J)=Mf/CFNRhow(Tx>)*A1) 

3200! 

an CALCULATE INSIDE AND OUTSIDE COEFFICIENTS 
3220! 

3230 Reu=FNRhHou(Tx)*Vu(J)*Di/FNMuu(Tx) 

3240 Cf=1, 

3250 Q=HMf<FMCpu(Tx)*w=(To(J)-Ti(J)) 

3260 Qp=0/(PI=*=Do*L) 

3270 IF (Mi=1 OR (Mi=2 ñND Нгуп<6)) AND J=0 THEN 
3280 Tfilm=Tsat 

3290 Kf=FNKu(Tfilm) 

3300 Rhof=FNRhow(Tfilm) 

3310 Hf g=FNHFg(Tsat)*1000 

3320 Muf =FNMuw(Tfilm) 

3330 Hnus.651*(Kf ^"3«Rhot ^2*Hfg*9.81/(Muf *Do«0p22^.32333 
3340 Tfilmc=-Tsat-Op/Hnu*.5 

3350 IF ABS((Tr1ilmc-Tf11m>)/Tf11mc>>.01 THEN 
3360 TfilmsTfi1lmc 
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3370 
3380 
3330 
3400 
3410 
3420 
3430 

RR" 
2440 
3450 
3460 
3470 
3480 
3490 
3500 
3910 
352U 
3530 


3650 
3660 
3670 
3680 
3630 
3700 
ЗА!" 
3720! 
3730! 
3740 


GOTO 3290 7 
END IF 
PRINT USING "ft0X.""Nusselt coefficient for first tube -» "",5D.D,"" (H/m^2) 


:Hnu 


END IF 

IF J=0 THEN 

IF Mi=1 AND Nrun=1 THEN Ho!=0, 

PRINT " Tube Vu Heat fiux Cond coef RI R2 


PRINT " z (m/S) (H/n^2) (4/m^2-K)" 
END IF - 
Мцш=Е ММцш С Тх )› 

Hi*FNKuCTx»/Di*Ci*Reu Ex*(FNPruCTx)2^.3333«Cf 
Dt-Qp/Hi*Do/Di 

Cfc=(Muu/(FNMuu(Tx+Dt))) .14 

IF BBS((OCF-Cfc)/Cfc)»5.01 THEN 

СЕе=(ФСЕ+С с)» ,5 

GUIU 34/U 

END IF 

Lmtd=(To(J)-Ti(J))/LOGC(Tsat-TitCJ))/(Tsat-TotJ))) 
Uo=Qp/Latad 

Ho(J)=1./¢01./U0-Do/(Di*Hi)-Rw) 

Rr=Uo/Ha( J) 

61-51-Но(.)) 

IF Nrun MOD S=1 THEN 

IF Мі-! QR <Mı=2 AND Nrun=31) THEN Ja=0 

IF Мі-2 AND S<Nrun AND Nrun<30 THEN Ja=Nrun-1 

IF Mi=2 AND 35<Nrun THEN Ja=Nrun-! 


IF Mi=1 OR (30<Nrun AND Nrun<36 AND Mi=1) OR Nrun<6 THEN 
R1=Ho(J)/H0(0) 

R2=S1/((J+1+Ja)*Ho(0)) 

ELSE 

Rt=HotJ)/Hol 

R2=Si/€CJ+1+Ja) *Hol) 

END ir 


PRINT RESULTS 
Poa USING "11X.DD, 4X,DD.D0,2X.2(D.5DE,2X),3(2.4D.2X)":;J+1+Ja.VutJ),0p,Ho 


СЕК „КЕ. Кт 


3750! 
3760! 
3770! 


IF Mi=2 AND Nrun<6 AND J=0 THEN 
Hol=Hol+Ho(0)/5 

END IF 

FOR K=0 TO 4 

IF K=J THEN S3(K)=S3(K)+R1 
IF K=J THEN S4(K)=S4(K)+R2 
NEXT K 

NEXT J 

IF run “0D S=0 THEN 

FOR K=0 704 

R3¢(K)=S3¢(K)/5 
RA(K)=S4CK) /5 

S3¢(K)=0. 

S4¢(K)=0. 

NEXT K 


‚ IF Mi=2 AND Nrun MOD 5=0 AND Nrun<>30 THEN Asi=Nrun#R4(4) «Но 1 


143 





2339: 
3960! 
3970 
3980 
3990 
4000 
4010 
4020 
4030 
4040 
4050 
4060 
4070 


РКІМІ нМехнНос KAILUS 


BRINT ет 

PRINT." Tube z R3 RAM 

FOR Js! TO 5 

PRINT USING "12X.DD.2€4X.Z.4D)":J*Ja, R3CJ- 12, RACJ- 1D 

OUTPUT 8F11e3:J*Ja, R3CJ- 12,R4CJ- 1) 

NEXT J 

EST "o" 

END IF 

гып MOD 5=0 AND Mi=2 AND Nrun<>30 AND Im=1 THEN 

PRINT USING "10X.""Set porous-tube flowmeter reading to '"'",3D.D.'" PERCENT 


vw... : Set 


4080 
4090! 
4100! 
4110! 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4130 
4200 


END IF 


IF Im=1 THEN 

BEEP 

INPUT "DO YOU HAVE ANY MORE DATA (1=YES,0=NO)2".Go_on 
er THEN Repeat 


IF Mi=2 AND Nrun<30 THEN Repeat 

IF Mi=1 AND Nrun<10 THEN Repeat 

END IF 

IF Im*1 THEN PRINT USING "!0X,DD,"" Data runs were stored in file "",10A": 


Nrun.NeudatasS 


4210 
4220 


ASSIGN 9»File! TO * 
ASSIGN eFile2 TO = 
ASSIGN 9*File3 TO + 
END 


THIS SUROUTINE CONVERTES THERMOCOUPLE VOLTAGE INTO TEMPERATURE 


SUB Tvsv (T) 
ОШЕА Сс) 
Sum=0. 

FOR I=0 TO 7 
Sum-Sum*C(CIO*T^I 
NEXT I 

T=Sum 

SUBEND 


THIS FUNCTION CALCULATES PRANDTL NUMBER OF WATER IN THE 


! RANGE 15 TO 4S DEG С 


DEF FNPru«(T) 

1-008961 7039766505-1*4(1.3749326E-2-T*(3.,968875E - 5-3. 45026E - 7*T»?)5) 
RETURN Y 

FNEND 


THIS FUNCTION CALCULATES THERMAL CONDUCTIVITY OF WATER 
IN THE RANGE OF 15 TO 105 DEG C 


DEF FNKu(T) 
Ү= ,5625894+Т7+ (2, 2964546Е-3-Т-+‹ 1 .509766Е-5-4,.0581652Е-3+*Т) ) 
RETURN Y : 


FNEND 
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5000! 
5010! 


5020! 
5040 
5050 
5060 
5070 
5030 
5090 
5100 
5110 


THIS FUNCTION CALCULATES SPECIFIC HEAT OF WATER 
IN THE RANGE 15 TD 45 DEG C 


DEF FNCpu«T? 
Y2€4,21120858-T*(2.26826E-3-T*(4.42261E-5*2.71428bE-7* T2552 +1000 


: RETURN Y 


FNEHD 


THIS FUNCTION CALCULATES DENSITY QF WATER IN THE 
RANGE 15 TO 105 DEG C 


DEF FNRhow( 7) 

Ко:399,52946+1+(.01269-Т» (5.482513Е-3-Т+1.234147Е-5)) 
RETURN Ко 

FNEND 


ШЕ асы гон APPLIES CORRECTIONS TO THERNSCOUPLE READING 


DEF FNTemp(T,D 

DIM A(14).B(14) 

ШРШ ЕЯО>535.2,573054,0.593101,0,57298.0.,56228.0.567354.0.569577/ 

ШАТА 0.553351.0.552008.0.566955,0.520998,0.522661.0.531008,0.560783.0.5524 


DATA 11.8744,8.63163.9.39412,3.570246.3.299436.3.36677.2.04507.7.459766 
ШІ /.4393225.7.9408,5.870/72,5.331556,6.13399,6.48585.,5Б.526224 

КЕАр А‹ +) ,В‹ +) 

SIS THEN 

ВН 0541) +.001=Т) 


FNEND 

THIS FUNCTION COMPUTES THE SPECIFIC VOLUME OF STEAM 

DEF FNYust(T) 
V=58.4525588-T+(1,51508776-T1+*(.01372746585-T*4,25366711E-5)) 
RETURN V 

FNEHD 


THIS FUNCTION CONVERTS THE VOLTAGE READING OF THE PRESSURE 


! TRANSIDUCER INTO PRESSURE IN MM HG 


DEF FNPvsv«(VU» 
Y=1.1302462+163.36413*V 
RETURN Y 

FNEND 


THIS FUNCTION CALCULATES THE SATURATION TEMPERATURE OF STEAM AS A FUNCTION 
RESSUR 


DEF FNTvsp«P) 
IF P«600 THEN 
T231.8776158*P*€.,235854929-P*(3.6613664E-4-P*2.41652372E -7)) 


ELSE 

[759 36562 *P« (07379467 -P* (3. 15662E -5-P «6 .27246E-9)) 
ND IF 

RETURN T 

FNEND 
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THIS FUNCTION COMPUTES THE VISCOSITY OF WATER 


DEF FNMw<T) 
Muse1.57609473E-3-T*(3.511938575E-5-T*(3.5835816E-7-1.3655861 15E -9« T2) 


- RETURN Mu 


FNEND 


t THIS FUNCTION COMPUTES THE LATENT HEAT OF VAPORIZATION 


DEF FHH£qCT) 
НҒа-2497.7389-Т%(2.2074%Т7%<1.7079Е-3-2.8593Е-6,Т)) 
RETURN Hfg 

FNEND 
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B. MODIFIED WILSON-PLOT METHOD 


10001 + nn EEE EEE TUE EE EE ЗЕ ЗЕ ЭЕ «ЗЕ Зе Є ЗЕ ЗЕ ЗЕ ЗЕ ЗЕ 


d FILE NAME: WILSON * 
1 ! ж 
1030! THIS PROGRAM COMPUTES THE SIEDER-TATE S 
1040! COEFFICIENT FOR FLOW IN TUBES м 


1050 Í ææ æ z x y æ e 36 -w- 36. 36 3€ 36 36. 3£ 36. 3 36 36. 36 2- 36 < 36. 38. 3€ 38. 36. 36 36 3 36. 3 36 36 3⁄ 36 3 3. -% 36 3 36 36 3⁄ 3 36 36 3 36 3 36 36 3£ 3 % Е 
ИЕШЕ СОМ /Сс, С‹?)› | 
ВОТ Yutl8) .Tc1 033) .Teo( 18) .Ts¢13) .Md¢ 18) .Taci8) 

1080 DATA 0.10086091 .25727 .94369, -767345 .8295. 78625535. 33 

1090 DATA -3247486589,6.97688b*11,-2.66192E*13,3.94078E*14 

1100 READ CC» 

PRINTER IS 701 

РЕЕР 


CLEAR 709 
INPUT "ENTER MONTH, DATE AND TIME (MM:D0:HH:MM:SSD",BS 
QUTPUT 709:"TD':Bs 
Di=.0141 
Ai=PI=Di*2/4 
Do=.915875 
2,305 
1200 Km=21.9 
1210 Ru=(Do-Di)*Do/(Km*(Do+Di)) 
1220 Rfis50. 
1230 Series: ! 
1240 OUTPUT 709;"TD" 
1250 ENTER 703;AS 
1260 PRINT USING "10X,""Month. date .and tıme: "",15A":AS 
1 ВЕЕР 
1280 INPUT "ENTER INITIAL GUESS FÜR SIEDER-TATE COEFFICIENT",Ci 


1300 ІНРОТ “ENTER EXPCNENT FOR REYNOLDS NUMBER" , Xn 
В RINISUSING "TOX.""Initial auess for Sieder-iate coefficient 
1320 PRINT USING "10X,""Exponent for the Reynoids number 


1340 INPUT “ENTER THE INPUT MODE (1=3054A.2=FILE)”.Im 

1350 IF Im=1 THEN 

1250 BEEP 

1370 INPUT “GIVE A MAME FOR THE DATA FILE".D fileS 

1380 CREATE BDAT D_files,5 

1330m ELSE 

1400 BEEP 

1410 INPUT “GIVE THE NAME OF THE DATA FILE".D fileS 

1420 PRINT USING "10X,""Foliouing analysis was performed for data stored ın fıl 

BAD Files - 

1430  BEEP 

1440 INPUT “ENTER THE NUMBER OF DATA RUNS STORED" .Nrun 

1450 END IF 

1460 ASSIGN 3File TO D fileS 

1470 BEEP 

1480 INPUT “GIVE A NAME FOR PLOTTING DATA FILE" ,DolotS 

1490 CREATE BDAT Dplots.5 

1500 ASSIGN sFi1lep TO DplotS 

15107 BEEP 

1520 INPUT ENTER ANALYSIS TYPE «TeHT,29UID»".Tt 

1530 IF It=1 THEN PRINT USING "10X,""Analysis type * Hi-METHOD''"" 
1540 IF [t=2 THEN PRINT USING "10X," "Analysis type = UL-METHOD""" 


шо со N ODO A Б ог) 
00050500000 


DIOS 
ESZ: DD An 


ии 
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1550 
1560 
1570 
1590 


К-0 

ES 

PRINT " Data Uu Tsat т То“ 
PRINT " е (m/S) (C) (C) cC» 


1590 Repeat:! 


1600! 
1610! 
1620! 
1630! 
1640! 
1650 
1660 
1670 
1680 
1650 


. 9 4 »^ 9 


2130 
2140 


RECORDS THERMOCOUPLE AND PRESSURE TRANSDUCER 
READINGS AUTOMATICALLY THROUGH THE HP 30548 
AUTOMATIC DATA ACQUISITION/CONTROL SYSTEM 


IF Im=1 THEN 

BEEP 

INPUT "ENTER FLONMETER READING (AS A PERCENT)" ,Fm 
OUTPUT 709:"AR AFO ALO" 

DUTPUT 722:"F1 R1 T1 Z} FLI” 

OUTPUT 709:"45 58" 


Бо АСПОЗ пи МОСОЧРЕ FOR WATFR INLET 


ШЕШКЕН” /22:1(0) 

OUTPUT 709:"AR AF3 RLS" 
DNI 22221 R1 711 Z1 FL1” 
FOR i-! TO 3 

OUTPUT 709:"AS SA” 


! READS THREE THERMOCOUPLES FOR WATER OUTLET 


ЕМШЕК 722:Т 1) 

NEXT I 

OUTPUT 209; AR AFIS ALIS" 
C УРТ КІ ТУТАТ РЕТ" 
DUTPUT 709; AS Sa” 


! READS PRESSURE TRANSDUCER 


ENTER 722:P_volts 
Pvap=FNPvsv(P_volts) 
Ts(K)=FNTvsp(Pvap) 
бит+0, 

ҒОК І-0 70 3 


! CONVERT VOLTAGE READINGS TO TEMPERATURE 


CALL Tvsv (TCI) 

ТЕ 150 THEN 
TeitK)=FNiemp(T(0),0) 

ES 

М=1+2 

APPLY ÍHERMOCOUPLE CORRECTIONS 
To=FNTemp(T(I),M) 

Sum=Sum+To 

END IF 

NEXT I 

Tco(K)=Sunm/3. 

ЕС ЗЕ 

ENTER ЗЕ 112: 715 СК) , Т 21 СК ) , ТсосКк) „Еп 
END IF 

Та(К) - <(Тс: (К) + Тсо<(К)?)».,5 
Md(K)=665.86*Fm/(100>650) 


! COMPUTE RATER-SIDE VELOCITY 


VwutK)=MdCK)/CFNRhou(TatK))*Aı) 

BEEP 

IF Im*1 THEN INPUT "ARE YOU TAKING MORE DATA «1-ҮЕ5.,0-М022",60 оп 
M=K+1 

PRINT USIHG "10X.DD.512X.,.DDD.DD)":M,VucK) ,Ts(K) ,‚Tcı(K) ,‚Tco(K) 
m. THEN OUTPUT &File:;Ts(K) ,Tcı(K) ,Tco(K) ‚Ffm 

K=K+ 


IF Im-1 THEN 
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2150 
2160 
2170 
2180 
2190 
2200 


2220! 


2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 


IF Go ons! THEN Repeat 
ELSE 

IF M<Nrun THEN Repeat 
END IF 

К-К- 1 

2-0 

21-0 


ДАЉИНА ItERATIGN TQ COMPUTE SIEDER-TATE COEFFICIENT 


Ssq=0 

5х-0. 

Sy=0. 

Sxs=0. 

Sxy=0. 

J=J+1 

BEIC ° 

PRINT " Iteration number = “;J 

D M OR Jj*1 THEN PRINT " X IUE 


FOR I=0 TO K 
O=Md(1>+«(Teo(T)-TeitI>)-4180., 
Lmtd=(Tco(I)-Tci(I))/LODG((Ts(I)-Tci(I))/(Ts(I)-Tco(I))) 
Un=0/(Lntd*PI=Do+L) 

Unr=1./Un 
Rei=FNRhout Tat I))+«VuCI)*=Di/FNMu(TatI?) 
Pru=FNPruw(Tati)) 

X=Rei `(-Xn)*Pru (-.3333) 

IF It=2 THEN 

Hir=Unr*D1/Do 

Ко: ЕНКъ(Та( ГТ) ) 

GOTO 2660 

END IF 

IF I=0 THEN 

Cf=1, 

Kw=FNKwtTa(I)) 
HisKw/Di*Ci*Rei^Xn*Pru^.2333*Cf 
Dt=Q/(PI«Di-L«Hi) 
Cfe=(FNMutTati))/FNMutTalId+Dt))".14 
IF ABS«KC(Cfc-C£»/CfFc)».01 THEN 
Cf=(Cf+Cfc)/2. 

GOTO 2500 

END IF 

Ho=1./(Unr-Rw-Rf 1*Do/Di-Do/(Di*Hi)) 
Нїг=1/Н1 

Оо=0 

ELSE 


CUM TE 1/HI 


Hir=(Unr-i/Ho«(Q/Qo) ~.3333-Ru-Rf 1*Do/Di)*Di/Do 
Dt=O0/(PI*Di*L)*Hir 

Cf =(FNMutTacI»»)/FNMutTa(i>+Dt>))7.14 
END IF 

И THEN 

Cf=1, 
Ніс=Кш#Сі/0 Веі ^Хп+Ргу”. 3333 + СҒ 
Dt=G/(PI*wDi*w*L*Hic) 
Cfc=(FHHu(Ta(I))/FNMudTa(I)+Dt)) ”.14 
IF ñRES((Cf-Cfc)/Cfc>.01 THEN 
(пар р- Сто)» 5 

GOTO 2680 
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2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 


END IF 

END IF 

Qp=O0/(PI*Do*L) 

X=X/Cf 

IF Jj=1 THEN 

OUTPUT @Filep;X,Hır 

Нігс-АзВжХ 

559=554+(Ніг-Нігс) ^2 

END IF 

ПНЕ THEN PRINT USING "10X.Z.7D,3X.2.7D,4X,DD.DD.2X.D.3DE.X.Z.53D 


К, г, ОЕ, бр, СЕ 


28401 


3060 
3070 
3080 
30930 
3100 
3110 
3120 
3130 
3140 
3150 
3160 


_files 


3170 
3180 


COMPUTE COEFFICIENTS FOR LEAST-SQUARES SCHEME 
Sx=Sx+X 

Sy=Sy+Hir 

Sxs=SxstX72 

Sxy=Sxy+Xx*Hir 

NEXT T 

МЕК + | 


* COMPUTE THE SLOPE QF THE LEAST-SQUARES LINE 


DEVELOPED FOR THE WILSON PLOT 

B=(N*Sxy-Sy*Sx)/(N#*Sxs-Sx*Sx) 

Cic=D:i/¢(B*Kw) 

PRINT USING "10X,""Intermediate value of Sieder-Tate coefficient = "",Z.4D 


IF ABS((Ci-Cic)/Cic)>.01 THEN 

Ci=Cic 

GOTO 2260 

END IF 

IF Jj=? THEN 

PRINT USING "10X.""Sieder-Tate coefficient 

А-(5у-Вж5х»-/М 

PRINT USING "10X.""Estimated foulıng factor 
PRINT " Least- squares iıne;" 

EMIT USING “13X.""Slope E SUB 

PRINT USING "13X,""Intercept = “MZ,SDE":A 
ELSE 

Jj=1 

GOTO 2260 

END IF 

PRINT USING "10X.'""Sum of squares = "“"",D.5DE":;:Ssq 
ASSIGN wFile TO * 

ASSIGN Filiep TQ = 

IF Im=1 THEN 


BEEP 
BRINT USING 210X, NOTE: "DD, =n Data runs were stored in-fule "",.8A";M,D 


END IF 
ОБ НС “10X.°° NOTE: “*,DD.°" X-Y pairs uere stored in file "",108";K* 


oa Сте 
МАП ЭШЕ,” 3 В 


1,Dplots 


3130 
3200 
3210 
3220 
3230 
3240! 
3250! 
3260! 
3270 
3280 
3290 


ВЕЕР 

INPUT "ARE YOU RUNNING ANOTHER SERIES (1=YES,0=-NO)?" бо оп 
IF Go_on=1 THEN Series 

END 

DEF FNRhow(T) 


THIS FUNCTION COMPUTES THE DENSITY OF WATER 
Ко-1006.35/24-Т»(.774489-Т“(2.262459Е-2-Т-«3.03304Е-4») 


RETURN Ro 
FNEND 
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3650 
3660 
3670 
3680 
3690! 
3700! 
3710! 
3720 
3730 
3740 
3750 
3760 
3770! 
3780! 
3790! 
3800 
3810 
3820 
3830 
3840! 
3650! 
3860! 
3570 
:3880 
3890 


DEF FNPvsv(V) 
THIS FUNCTION CONVERTS THE PRESSURE TRANSDUCER 


! READING FROM VOLTS TO PRESSURE IN MM HG 


Y=1.1103462+163.36413*V 
RETURN Y 

FNEND 

DEF FNTvsp(P) 


THIS FUNCTION COMPUTES THE SATURATION TEMPERATURE 


! CORRESPONDING TO PRESSURE IN MM HG 


IF P«600 THEN i 

А а F-(9:6613664E- 4-P»2.41652372E- 7) 
Т-59,36552%Ре< ,0737945/-Р»(3.15552Е-5-2-5.27246Е-9); 

END IF 

RETURN Т 

FNEND 

DEF FNTemp(1,D 


THiS FUNCTION APPLIES THERMOCOUPLE CORRECTIONS 


DIM ACS).BCS) 

DATA 0.640533.0.573054,0.593101.0.57298.0.567284.0.569577 
DATA 11.9744,8.63163,9.398412.8.570246,8.299435.8.36677 
READ ас) ВС) 

T=T-(ACTI)-BCI)*.001*T) 

RETURN T 

FNEND 

SUB Tvsv(T) 

ОШМУ Сс/ C¿e7) 

Sum=0. 

FOR I=0 TO 7 

Sum=Sum+C(I)*T`I 

NEXT I 

T=Sum 

SUBEND 


THIS FUNCTION COMPUTES PRANDTL NUMBER FOR WATER 


DEF FNPru(T) 
Prs10'(1.09976605-T*C.013749326-1*(3.968875E-5*43.45025E-7* 120) 
RETURN Pr 

FNEND 

DEF FNMu«T) 


THIS FUNCTION COMPUTES THE VISCOSITY OF WATER 
Mus1.5087546575E-3-T*(3.025732489E-5- T* (2.626439826E-7-T*8.18601927E- 100) 
RETURN Mu 

FNEND 

DEF FHKu( T) 

THIS FUNCTION COMPUTES THE THERMAL CONDUCTIVITY OF WATER 

Kwa= .572183504477+1.52770121209E-3+T 


RETURN Kua 
FNEND 
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Ge 


1000! 
1010 


1020: 


1030 
1040 


— <> осо мол љом — 
O D O <> <> © < > <> < <> 


— — — — — — —— — — 
№“ М) — — mu о — — 
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PLOTTING PROGRAM 


FILE NAME: PLOT 

PRINTER IS 705 

БЕЕР 

qudd “ENTER MINIMUM AND MAXIMUM X-VALUES” .Xmin.Xmax 

Е "ENTER MINIMUM AND MAXIMUM Y-VALUES’. Ymin. {max 

ш "ENTER STEP SITZE RüReX-OXIS'sXs tep 

N "ENTER STEP SLZE.FDR Y-AXIS”.Ystep 
E - 

PRIMT "IN:;SP!:IP 2390.1800. 9300. 6800;°° 

PRINT “SC 0. 100. ОССЕ 2,0: 

STAT TOUS Satan Alta ter 

Sfy=100/(Ymax-Ynin) 

PRINT “PU 9.0 PD” 

FOR Xa=Xmin TQ Xmax STEP Xstep 

X=(Xa-Xmin) «St x 

КИТА РА: Х "0: ХТ:"“ 

МЕХТ Ха 

PRINT "PA 100.0:PU:" 

PRINT “PU PA 0,9 PD” 

FOR Ya=Ymin TO Ymax STEP Ystep 

Y=(Ya-Ymind*«Sfy 

КЕЧЕ РА 0. :Y. YT" 

NEXT Ya 

PRINT “РА 0,100 TL OQ 2" 

FOR Xa=Xmin TO Xmax STEP Xstep 

Х- СХа-Хпіп)ж5Ғх 

ЕКЕЛШ X." .160: XT" 

NEXT Xa 

PRINT "PA 100,100 PU PA 100.0 PD" 

FOR Ya=Ymin TG Ymax STEP Ystep 

Y=(Ya-Ymind#Sfy 

ШИН ro РА 10066" .%. УТ" 

NEXT Ya 

PRINT "PA 100,100 PU" 

poe en 0.-2 SR 1.5.2" 

FOR Xa=Xmin TO Xmax STEP Xstep 

X=(Xa-Xmin) «Sfx 

ЕСО РАУ tt. 0." 

ИСР -2.-1;LB':Xa;" 

NEXT Xa 

PRINT “PU РА 0.0" 

FOR Ya=Ymin TÜ Ymax STEP Ystep 

Y=(Ya-Ymin)*Sfy 

PRINT “PA O0, "ү, Sd 

ПРОПИСА =аА = 25:18 :Ха;"" 

NEXT Ya 


INPUT "ENTER X-LABEL'".,XlabelS 
ТУР ENTER Y-LABEL”.Ylabels 
PRINT "SR T.S.2:PU РА 50.-10 CP'":-LENCX1abelS)/2:;"0; ae 


ГАШ “РА -11‚50 CP 0. ;-LEN(YlabelS)/2*5/6:' ИОС ЕВИ 
PRI CP 0,0 DI" 
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1560 
1570 
1580 
1590 
1600 
1610 


1620 


1§30 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1210 
1720 
1755 
1740 
1750 
1760 
1220 
1780 
1790 
1800 
1810 
1820 
1830 
1940 
1850 
1860 
1870 
1880 
1830 
1900 
1910 
1920 
19 0 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2190 
2140 
2150 


Repeat:! 


ВЕЕР 

ШЕШІ ПП ҮШ) НВЫТ І0 РЕОТ DATA FROM A FILE (1+YES.0=NO) 2?” .Ok 
IF Ok=1 THEN 

ВЕЕР 

DIT ЗЕМТЕК THE NAME OF THE DATA FILE".D.fileS 

ASSIGN “File TO D_files 


BEER 

INPUT "ENTER THE BEGINNING RUN NUMBER” ,Md 

ВЕЕР 

INPUT “ENTER THE NUMBER OF X-Y PAIRS STORED" .Npairs 
ВЕЕР 


INPUT "SELECT A SYMBOL FOR THE PLOTTER (1=*,2=+,3=c0,4=0,5= )",5y 
PRINT "PU DI" 
IF Sy=1 THEM PRINT "SM«' 
IF 5у=2 THEN PRINT "SM" 
IF Sy=3 THEN PRINT "5Мс“ 
EE EN PRIMI "SMs' 
IF Sy-5 THEM PRINT "SM^'" 
BEEP 
INPUT "SELECT MODE (1=HN/H1.2=HMNC(avg)/H1)” , Jm 
IF Md>Npairs THEN 
FOR f=! TO (Md-1) 
IF Jm=1 THEN ENTER #File:Xa.Ya,Yy 
IF Jm=2 THEN ENTER %Fıle;Xa,Yy,Ya 
NEXT I 
END IF 
РОК 121 10 Npairs 
IF Jm=! THEN ENTER &Fi1le:Xa,Ya,.Yy 
IF Jm=2 THEN ENTER ®File:Xa,Yy.Ya 
X=(Xa-Xmin)<«Sfx 
Y=(Ya-Ymin)d#Sfy 
PRINT "PRU. X.Y, 
NEXT I 
BEEP 
ASSIGN File TO = 
INPUT “DO YOU HAVE MORE DATA TO BE PLOTTED (1=YES,0=NQ) 2?" бооп 
IF Go_on=1 THEN Repeat 
END IF 
ВЕЕР 
INPUT “DO YOU LIKE TO PLOT THE NUSSELT RELATION (1=YES,0=N0)?",Go_on 
PRINT “PW SM" 
IF Go_on=! THEN 
FOR Xa=Xmin TO Xmax STEP Xstep/50 
X=(Xa-Xmin)*Sfx 
IF Jm=1 AND Xa>Xmin THEN Ya-Xa^.75-(Xa-12^.75 
ТЕ Jm=2 AND Xa»Xmin THEN Ya-Xa (-.25) 
IF Xa=Xmin THEN Ya=1 
Y=(Ya-Ymın)*Sfy 
Pet  PA’,X,.Y, PD" 
NEXT Xa 
BEES 
PRINT “PU” 
INPUT “MOVE THE PEN TO LABEL THE NUSSELT LINE".Ok 
PRINT "LBNusselt" 
END IF 
BEEP 
INPUT “DO YOU LIKE TO PLOT KERN RELATIONSHIP?”. Yes 
IF Yes=1 THEN 
FOR Xa=Xmin TO Xmax STEP Xstep/20 
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2440 
2450 


2470 
2480 


Та-Ха”(-1/Б» 
X=(Xa-Xmin) «Sfx 
Y=(Ya-Ymin)*Sfy 
PRINT "PA",X.Y,"PD" 


NEXT Xa 
РАТЕ "РО" 
ВЕЕР 


INPUT “MOVE THE PEN TO LABEL KERN RELATIONSHIP", Ok 
PRINT "LBKern:PU" 

END IF 

PRINT "PU Pñ 0.0" 

BEEP 

INPUT “OK TO PLOT CURVE FOR NEW DESIGN (1=0K,9=N0)?" .Dk 
IF Ov =1 THEN 

FOR Xa=Xmin TO Xmax STEP Xstep/2 

IF Xa=1 THEN Ya=1 

IF Xa>i THEN Ya=((Xa*5)” .75-(Xa*5-1)".75+5)/6 
X=(Xa-Xmin)*Sfx 

Y=(Ya-Ymin)«Sfy 

РЕШИ "PA”.X.Y."PD" 

МЕХТ Ха 

PRINT "Py" 

END IF 

BEEP 

INPUT "DO YOU LIKE TO PLOT EISSENBERG RELATION?",Go on 
IF Go. on*1 THEN 

FOR Xa=Xmin TO Xmax STEP Xstep/10 

(Не Бо” а2яХа (-,25) 

X=(Xa-Xmin)*Sfx 

Y=(Ya-Ymin)*Sfy 

PRINT "PA",X,.Y ."PD" 

NEXT Xa 

PRINT "PU" 

BEEP 

INPUT "MOVE THE PEN TO LABEL THE EISSENBERG LINE" .Ok 
PRINT "LBEissenberqg;PU" 

END IF 

PRINT “PU” 

БЕЕР 

INPUT "DO YOU LIKE TO PLOT THE EXPTL CURVE (15T,02NO0»?",Go on 
Eo nl THEN 


INPUT "ENTER THE СХРОМЕМТ” ‚Ех 

FOR Xa=Xmin TO Xmax STEP Xstep/10 

Уа=Ха`(-Ех) 

X=(Xa-Xmin) «Sfx 

Y=(Ya-Ymin)*Sfy 

PRINT "PA".X,Y,°FD". 

NEXT Xa 

PRINT “PU” 

BEEP 

INPUT "MOVE THE PEN TO LABEL THE EXPTL CURVE”, Ok 
PRINT "LBHNCavg) /H1=N" 

PRINT PR 1.1” 

Peer EB": -Ex;"" 

END IF 

BEEP 

INPUT "DO YOU LIKE TQ DRAW A STRAIGHT LINE?" ,бо_оп 
son! THEN 


INPUT "ENTER THE SLOPE",SI 
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292 

2930 
2940 
2950 
2950 
2970 


BEEP 

INPUT "ENTER THE INTERCEPT".Ac 
FOR Xa=Xmin TO Xmax STEP (Xmax-Xmin? 
Үа=Ас+5 1 +Ха 
Y=(Ya-Ymin)*Sfy 
Х-(Ха-Хпіп)5Ғх 

IF Y«0 THEN 

Хат= Утап-ас)/51 
X=(Xam~-Xmin) «Sfx 
Y=0 

END IF 

IF Y>190 THEN 
Xam=(Ymax-Ac>/Sl 
X=(Xam-Xmın)*Stx 
Y=100 

END IF 
ERINTT SPA" .X.Y."PD" 
NEXT Xa 

END IF 

PRINT “PU SPO" 

END 
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